GO 


AFAPL-TR-76-57 


* 


HEAT  ADDITION  TO  A SUPERSONIC  FLOW 


RAMJET  APPLICATIONS  BRANCH 
RAMJET  ENGINE  DIVISION 


OCTOBER  1976 


D D C 

-?.( — >rpnp  nrp 


NOV  23  1976 

ulkiaoTTE 


TECHNICAL  REPORT  AFAPL-TR-76-57 

FINAL  REPORT  FOR  PERIOD  1 SEPTEMBER  1972  to  1 JANUARY  1974 


Approved  for  publi  relea-.,  ' stribution  unlimited 


AIR  FORCE  AERO  PROPULSION  LABORATORY 
AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a definitely  related  Government  procurement  operation, 
the  United  States  Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever;  and  the  fact  that  the  Government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be 
regarded  by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture, 
use,  or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  has  been  reviewed  by  the  Information  Office,  (ASD/OIP)  and  is  re- 
leasable to  the  National  Information  Service  ( NT IS).  At  NTIS,  it  will  be 
available  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 

dP'/?.  C 7]utstA-c*\  jf*  ■ 

DR.  WILLIAM  C.  BURSON,  Jr.,  GS-13 
Aerospace  Engineer 
Ramjet  Applications  Branch 
Ramjet  Engine  Division 


FOR  THE  COMMANDER 


Strategic  Missile  Propulsion 
Ramjet  Applications  Branch 
Ramjet  Engine  Division 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by 
security  considerations,  contractual  obligations,  or  notice  on  a specific  document. 


SECURITY  CLASSIFICATION  of  This  PAGE  'When  Data  Entered, 


READ  INSTRUCTIONS 


REPORT  DOCUMENTATION  PAGE 


BEFORE  COMPLETING  FORM 


2 GOVT  ACCESSION 


AFAPL-TR- 76-57 


Heat  Addition  to  a Supersonic  Flow# 


D*\  William  C./Burson,  Jr 


9 PE  RFOPMINO  ORGANI  7 AT|QN  NAME  AND  ADDRESS 


AFAPL/RJA 

Wright-Patterson  AFB,  Ohio  45433 


2h403 


" CONTROLLING  OFFICE  MIME  UNO  ADDRESS 


Air  Force  Aero  Propulsion  Laboratory  ( / / j 

Ramjet  Engine  Division,  AFAPL/RJA 
Hright-Patterson  AFB,  OH  45433 

U MONI  TQPING  AGENCY  NAME  ft  ADDRESS'!/  different  from  Controlling  Otl  e 


Approved  for  public  release;  distribution  unlimited 


17  DISTRIBUTION  STATEMENT  < of  the  abstract  entered  In  Block  2".  II  different  from  Report) 


19  KEY  WORDS  ^Continue  on 

Supersonic  Flow 
Jet  Propulsion 
Plasma  Flow 
Gas  Dynamics 


■erse  aide  if  necessary  end  identify  by  block  number. 


reverse  side  If  necessary  and  identify  by  block  number) 


O^T  RACT  / Continue 

Advanced  propulsion  engines  of  the  future,  which  may  fly  at  Mach  numbers 
greater  than  five,  will  depend  upon  the  supersonic  heat  addition  process  for 
efficient  operation.  In  supersonic  heat  addition,  the  gas  flow  remains  super- 
sonic before,  during,  and  after  the  heating  process.  The  advantages  of 
supersonic  heat  addition  compared  to  subsonic  heat  addition  are  a lower  static 
temperature  and  a lower  static  pressure  for  a given  enerqy  input. 

Two  supersonic  heat  addition  modes  are  available.  One  is  the  subsonic 
diffusion  flame  propagating  into  a supersonic  flow  where  the  normal  or  — *'  ^ 


FORM 
1 JAN  73 


EDITION  OF  1 NOV  65  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  This  »AGE  't*"hen  Data  Entered] 


Unclassi fied 


SE  -uRlTY  CLASSIFICATION  OF  ThIS  PAGE  »'hon  Data  En  farad) 


'perpendicular  component  of  the  flow  velocity  with  respect  to  the  flame  front 
is  subsonic.  The  disadvantage  of  the  subsonic  diffusion  flame  is  the  relatively 
long  combustor  length  resulting  from  the  subsonic  normal  component  of  the 
flow  velocity.  The  other  supersonic  heat  addition  mode  is  a weak  detonation 
wave  where  the  thermal  wave  or  flame  front  is  supersonic  relative  to  the  heated 
and  unheated  gas.  The  purpose  of  this  effort  was  to  experimentally  investigate 

the  weak  detonation  wave  by  simulating  it  in  a gas  flow  and  a water  flow.x^ 

The  gas  flow  weak  detonation  wave  was  produced  by  radio  frequency  (R.F.)  heating 
of  a supersonic  Argon  flow  using  an  R.F.  coil  wrapped  around  a Boron  Nitride 
duct.  The  water  flow  weak  detonation  wave  was  simulated  by  mass  addition  to  a 
supersonic  water  stream. 

A quasi-one-dimensional  theory  based  on  average  properties  was  developed 
which  can  analyze  axisymmetric  duct  flows  with  wall  friction,  varying  area, 
heat  addition,  and  mass  exchange  occurring  simultaneously.  The  theory  was 
used  for  comparison  with  the  experimental  results. 


Unclassified 


AFAPL-TR-76-57 


FOREWORD 

This  report  contains  the  results  of  an  investigation  of  heat 
addition  to  a supersonic  flow.  The  work  was  performed  in  the  Ramjet 
Engine  Division  of  the  Air  Force  Aero  Propulsion  Laboratory,  Air 
Force  Systems  Command,  Wright-Patterson  AFB,  Ohio,  under  Project  3012, 
Task  14,  and  Work  Unit  03.  The  work  effort  was  conducted  by  William  C. 
Burson,  Jr.,  during  the  period  September  1972  to  January  1974.  Mr. 
Stanley  R.  Double  of  AFAPL/TFP  served  as  technician  for  the  effort. 

The  author  wishes  to  express  here  his  gratitude  to  Professor  R. 
Edse  for  the  guidance  and  encouragement  received  as  his  advisor  in 
the  doctoral  program  at  The  Ohio  State  University.  The  author  wishes 
also  to  give  recognition  to  the  Air  Force  Institute  of  Technology, 
particularly  to  Dr.  Andrew  Shine  and  Dr.  Harold  Wright  of  the 
Aerospace-Mechanical  Engineering  Department,  who  permitted  the  use 
of  their  water  table  test  facility;  to  Mr.  Stanley  Double,  the  tech- 
nician at  AFAPL  for  his  able  assistance  in  preparing  and  performing 
| many  of  the  experiments;  to  Mr.  Kenneth  Johnston  of  the  Air  Force 

Computer  Center  for  his  assistance  in  preparing  and  debugging  the 
computer  programs  used  in  this  effort;  and  to  Miss  Cynthia  Lines  for 
the  exceptionally  well  handled  typing  task. 

iii 


I 


AFAPL-TR-76-57 


TABLE  OF  CONTENTS 


SECTION 

PAGE 

I 

INTRODUCTION 

1 

II 

BACKGROUND 

4 

III 

THEORETICAL  ANALYSIS 

12 

Flow  Equations 

12 

Water  Table  Theory 

41 

Electric  Probe  Theory 

44 

R.F.  Energy  Addition 

60 

IV 

EXPERIMENTAL  EQUIPMENT  AND  PROCEDURE 

64 

Gas  Flow  Apparatus 

64 

Flow  Measurement  Probes 

69 

Water  Table  Apparatus 

78 

V 

EXPERIMENTS 

83 

Water  Table  Experiments 

83 

Gas  Flow  Experiments 

90 

VI 

SUMMARY  AND  CONCLUSIONS 

129 

Experiments 

129 

Radial  and  Axial  Variations 

130 

Average  Properties  for  Gas  Experiments 

133 

Conclusions 

135 

LIST 

OF  REFERENCES 

141 

PRECEDING  PAGE  ELANK-NOT.  FILLED 

" II  ""  ■« 


AFAPL-TR-76-57 


15  Schematic  of  the  Langmuir  Probe  Characteristic  Curve  ....  46 

16  Voltage-Current  Characteristic  of  the  Double  Probe  48 

17  General  Potential  Diagram  for  the  DPM 48 

18  Comparison  of  Experimental  and  Theoretical  DPM  Traces  ...  50 

19  Domain  of  Validity  for  Probe  Theories  With  Data  Points 


Shown 52 


AFAPL-TR-76-57 


LIST  OF  ILLUSTRATIONS  (CONTb) 

Figure  Page 

21  Test  Apparatus 66 

22  Photograph  of  Test  Apparatus  With  Power  Input  to  the  Gas  . . 67 

23  Photograph  of  Apparatus  for  Supersonic  Gas  Flow  Energy 

Addition  Tests  68 

24  Pressure  Probes 70 

25  Static  and  Total  Pressure  Probes  Along  Side  of  a 

Boron  Nitride  Duct 70 

26  Pressure  Measurement  Schematic  72 

27  Total  Temperature  Probe  Schematic 74 

28  Double  Probe 7 5 

29  Double  Probe  Schematic  76 

30  Schematic  of  Double  Probe  Circuit 77 

31  Schematic  of  Collector 79 

32  Collector  As  Seen  Inside  the  Vacuum  Chamber.  To  the  Right 
Is  the  Port  In  Which  the  Collector  Is  Placed  During  Gas 

Flow  Tests 79 

33  Water  Table  Schematic 80 

34  Water  Table  with  Nozzle  Direct-Connect  Duct  and  Porous 

Injection  Plate 81 

35  Water  Table  Test  with  Free  Inlet  Duct 81 

36  Water  Table  Results 84 

37  Side  View  of  Duct  Directly  Connected  to  Nozzle 85 

38  Mass  Addition  or  Heat  Addition  Greater  Than  in  Figure  37  . . 85 

39  Side  View  of  Free  Inlet  Duct  Without  and  With  Mass 

Addition  (Heat)  Addition  86 

40  Side  View  of  Free  Inlet  Duct.  Mass  Addition  Greater  Than 

in  Figure  39 87 

41  Side  View  of  Free  Inlet  Duct.  Mass  Addition  Greater  Than 

in  Figure  40 87 

vii 


AFAPL-TR-76-57 


LIST  OF  ILLUSTRATIONS  (CONTb) 

Figure  Page 

42  Side  View.  No  Mass  Addition  and  Ink  in  Water 88 

43  No  Mass  Addition.  Length  of  Duct  Extended  Illustrating 

the  Effect  of  Friction 88 

44  Static  Press,  vs.  Axial  Position  89 

45  Varying  Area  Duct  for  Heat  Addition 91 

46  Pressure  for  Wall  Taps  & Static  Pressure  Probe 95 

47  Pressure  Measured  by  Wall  Taps 96 

48  Pressure  Measured  by  Wall  Taps 58 

49  Probes  In  Flow  At  Duct  Exit 98 

50  Mach  Number  vs.  Radial  Position 101 

51  Mach  Number  vs.  Radial  Position 102 

52  Average  Mach  No.  vs.  Axial  Position 103 

53  Average  Press.  Ratio  vs.  Axial  Position 104 

54  Static  Press.  Ratio  vs.  Radial  Position 105 

55  Total  Press.  Ratio  vs.  Radial  Position  106 

56  Average  Total  Press.  Ratio  vs.  Axial  Position 107 

57  Static  Temp.  Ratio  vs.  Radial  Position  108 

58  Static  Temp.  Ratio  vs.  Radial  Position  109 

59  Total  Temp.  Ratio  vs.  Radial  Position 110 

60  Total  Temp.  Ratio  vs.  Radial  Position Ill 

61  Velocity  Ratio  vs.  Radial  Position  112 

62  Density  Ratio  vs.  Radial  Position 113 

63  Watts  Per  Sq.  Inch  vs.  Radial  Position 114 

64  Mass  Flow  Ratio  vs.  Radial  Position 115 

65  Total  Pressure  vs.  Total  Temperature  Ratio  for  Constant 

Area  and  Diverging  Ducts 136 


viii 


AFAPL-TR-76-57 


LIST  OF  ILLUSTRATIONS  (CONCL'D) 

Figure  Page 

66  Mass  Flow  Ratio  Caused  by  Heat  Addition  vs.  Total 

Temperature  Ratio  for  Constant  Area  & Diverging  Ducts  ....  137 

LIST  OF  TABLES 

Table 

1 T°  - Measured  By  Probe 93 

2-4  Comparison  of  Experiment  and  Theory 

Constant  Area  Duct 116 

5-7  Comparison  of  Experiment  and  Theory 

Area  Ratio  of  1 .60 119 

8-10  Comparison  of  Experiment  and  Theory 

Area  Ratio  of  2. 0 122 

11  Double  Probe  Data  For  Three  Total  Temperature 

Ratios  as  a Function  of  Radial  Position. Constant 
Area  Duct 


128 


AFAPL-TR-76-57 


I 


i 

LIST  OF  SYMBOLS  (CONCL ’D) 


r 

i 


RO 

S 

T 

T 

e 

T. 


1 


u 


V 


v 

X 


y 


$ 

p 

Di 

Y 

8 


a 

XD 

\ 

C 

P 

v 


Duct  Radius 
Entropy 

Static  Temperature 


Electron  Temperature 
Ion  Temperature 
Total  Temperature 
Total  Temperature 
Velocity 

Velocity,  Voltage 
Specific  Volume 
Distance  or  Displacement 

Mass  Injection  Velocity  Component  Parallel  to  Main  Flow 
Divided  by  the  Main  Flow  Velocity 


Compressor  Ratio 
4f  + 4m  (l-y)n 
Den  si ty 

rp/XD 

Ratio  of  Specific  Heats 

Divergence  or  Convergence  Angle 

Degree  of  Ionization 

Debye  Length 

Mean  Free  Path 

Mobi 1 i ty 

Collision  Frequency 
Power  Source  Frequency 
Permi tti vi ty 


x 


AFAPL-TR-76-57 


LIST  OF  SYMBOLS  (CONT ’D) 

k Boltzman's  Constant 

K.E.  Kinetic  Energy 

1 Length  of  Tungsten  Wire 

L Duct  Length 

LN  n 

Ln  Natural  Logarithm 

m 2 tane 

nr  Ion  Mass 

m Electron  Mass 

e 

m Neutral  Particle  Mass 

n 

m Mass  Flow  Rate 

M Mach  Number 

n Variable  Parameter 

n Ion  Particle  Density 

o J 

n Electron  Particle  Density 

e 

n Metastable  Number  Density 

m J 

N Variable  Parameter,  Power  Setting 

P Static  Pressure,  Power  Density 

Pow  Power 

P°  Total  Pressure 

PH  f/tane 

Q Specific  Energy  Input  to  the  Gas 

q Water  Mass  Injection  Term 

r Recovery  Factor 

rp  Probe  Radius 

R Radius 

xi 


d 


AFAPL-TR-76-57 


LIST  OF  SYMBOLS 

A Area 

c.j  Ion  Average  Velocity 

Cp  Specific  Heat  at  Constant  Pressure 

Cv  Specific  Heat  at  Constant  Volume 

D Diameter 

e Electronic  Charge 

E Effective  Applied  Electric  Field 

a 

Em  Energy  Difference  Between  Metastable  and  Ionization  Level 

EV  Recombination  Energy 

f Friction  Factor 

F Froude  Number 

g Gravity  Acceleration 

g Statistical  Weight  of  Metastable  State 
m 

g-  Statistical  Weight  of  Ion  Ground  State 
gg  Statistical  Weight  of  an  Electron 

G Ratio  of  Specific  Heats 

h Water  Height,  Enthalpy,  Planck's  Constant 

i Electric  Current 

I Ion  Saturation  Current 

I Axial  Duct  Position 

J Current  Density 


AFAPL-TR-76-57 


SECTION  I J 

INTRODUCTION  j 

Advanced  propulsion  engines  of  the  future,  which  may  fly  at  Mach  j 

numbers  greater  than  five,  will  depend  upon  the  supersonic  heat  addi- 
tion process  for  efficient  operation.  In  supersonic  heat  addition, 
the  gas  flow  remains  supersonic  before,  during,  and  after  the  heating 
process.  The  advantages  of  supersonic  heat  addition  compared  to  sub- 
sonic heat  addition  are  a lower  static  temperature  and  a lower  static 
pressure  for  a given  energy  input.  The  lower  pressure  and  temperature 
provide  a less  severe  structural  requirement  for  the  vehicle;  also, 
the  lower  temperature  results  in  less  dissociation  of  the  gas.  Disso- 
ciation is  an  energy  loss  if  recombination  does  not  occur  within  the 
engine. 

Two  supersonic  heat  addition  modes  are  available.  One  is  the 
subsonic  diffusion  flame  propagating  into  a supersonic  flow  where  the 
normal  or  perpendicular  component  of  the  flow  velocity  with  respect  to 
the  flame  front  is  subsonic.  The  disadvantage  of  the  subsonic  diffusion 
flame  is  the  relatively  long  combustor  length  resulting  from  the  sub- 
sonic normal  component  of  the  flow  velocity.  The  other  supersonic  heat 
addition  mode  is  a weak  detonation  wave  where  the  thermal  wave  or  flame 
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front  is  supersonic  relative  to  the  heated  and  unheated  gas.  Supersonic 
flame  speeds  are  virtually  impossible  to  obtain  by  thermal  ignition  but 
may  be  possible  by  photochemical  ignition.  The  advantage  of  the  weak 
detonation  wave  is  the  short  combustor  length  that  it  would  require. 

The  purpose  of  this  effort  was  to  experimentally  investigate  the 
weak  detonation  wave  by  simulating  it  in  a gas  flow  and  a water  flow. 

The  gas  flow  weak  detonation  wave  was  produced  by  radio  frequency  (R.F.) 
heating  of  a supersonic  Argon  flow  using  an  R.F.  coil  wrapped  around  a 
Boron  Nitride  duct.  The  coil  length  represents  the  wave  region  since 
the  heating  is  by  electrons  which  is  quite  rapid.  Comparison  of  the 
experimental  measurements  was  made  with  a theoretical  analysis.  The 
water  flow  weak  detonation  wave  was  simulated  by  mass  addition  to  a 
supersonic  water  stream.  The  water  flow  simulation  was  performed  to 
determine  if  it  could  qualitatively  predict  what  occurs  in  gas  flow 
heat  addition. 

The  water  flow  experiments  were  performed  on  a water  table  using  a 
two-dimensional  nozzle  to  simulate  a Mach  2 exit  flow.  The  "combustor" 
at  the  nozzle  exit  consisted  of  two  plexiglas  walls,  one  on  each  side 
of  a porous  injection  plate  that  allowed  mass  injection  to  the  main 
flow.  The  results  of  the  water  table  experiments  did  predict  the 
phenomena  seen  in  the  gas  flow  experiments. 

Supersonic  heat  addition  to  supersonic  Argon  gas  flow  was  achieved 
in  constant  area  and  varying  area  ducts.  The  low  Reynolds  number 
(35,000)  of  the  flow  caused  some  difficulty  in  obtaining  good  pressure 
measurements  and  permitted  considerable  upstream  influence  of  the 
pressure  rise  generated  by  heat  addition.  A pressure  rise  and  decrease 
within  the  duct  was  observed  during  heat  addition.  This  pressure 
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A 


variation  within  the  duct  could  be  caused  by  two  phenomena.  First,  a 
separation  region  and  reattachment  could  occur;  and,  second,  due  to  the 
non-uniform  radial  energy  addition,  the  rapid  heating  near  the  wall 
within  the  coil  region  causes  a rapid  expansion  next  to  the  wall, 
forcing  more  mass  flow  towards  the  center  which  effectively  decreases 
the  flow  area.  Downstream  of  the  coil,  radial  heat  conduction  occurs, 
thus  relieving  the  pressure.  Measurements  were  made  with  a static  pres- 
sure probe,  a total  pressure  probe,  a total  temperature  probe,  and  an 
electric  double  probe  for  determination  of  the  ionization  level.  The 
degree  of  ionization  had  to  be  known  in  order  to  make  a correction  to 
the  total  temperature  measurement. 

A quasi-one-dimensional  theory  based  on  average  properties  was 
developed  which  can  analyze  axisymmetric  duct  flows  with  wall  friction, 
varying  area,  heat  addition,  and  mass  exchange  occurring  simultaneously. 
The  theory  was  used  for  comparison  with  the  experimental  results. 
Satisfactory  agreement  between  theory  and  experiment  was  obtained  if 
the  average  of  the  measured  temperature  was  obtained  by  averaging  with 
a radial  mass  flow  variation  rather  than  just  averaging  with  the  radius. 
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SECTION  II 
BACKGROUND 

Heat  addition  to  a supersonic  flow  is  of  great  interest  in  the 
design  of  advanced  propulsion  systems.  A thermodynamic  analysis,  based 
on  the  Brayton  Cycle,  shows  that  for  the  maximum  exit  Mach  number  or 
thrust  of  a jet  engine  the  compressor  pressure  ratio  (:)  approaches  one 
when  the  flight  Mach  number  (Ma)  is  3.87  and  the  total  temperature  ratio 
(total  temperature  measured  after  heat  input  divided  by  the  ambien*  air 
temperature,  T2°/Ta)  is  four.  A total  temperature  ratio  equal  to  four  is 
based  on  limitinq  the  uncooled  turbine  blade  temperatures  to  2000  R for 
an  ambient  temperature  of  500  R.  Of  course,  cooled  turbine  blades  will 
allow  higher  flight  Mach  numbers  than  3.87  for  : * 1.0  at  the  maximum 
thrust  condition.  The  optimum  compressor  ratio  for  the  above  is: 


A compressor  ratio  approaching  one  indicates  that  a ramjet  instead  of  a 
turbojet  engine  s ""11  be  used.  The  ramjet  engine  consists  of  an  inlet, 
diffuser,  burner,  and  exhaust  nozzle.  In  the  Mach  3 to  5 region,  satis- 
factory operation  of  the  ramjet  is  possible  by  means  of  supersonic 
diffusion  throuqh  several  oblique  shocks  which  are  followed  by  a normal 
shock  and  subsonic  diffusion  so  that  energy  is  addea  by  subsonic  combus- 
tion. The  overall  efficiency  of  the  ramjet  at  Mach  numbers  of  5 and 
greater  depends  very  much  on  the  mode  of  heat  addition.  It  appears  that 
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the  most  advantageous  process  involves  diffusion  to  a lower  Mach  number 
which  is  still  supersonic  while  the  energy  addition  or  combustion  occurs 
in  this  supersonic  air  stream. 

The  advantages  of  supersonic  heat  addition  or  combustion  compared  to 
the  subsonic  case  are  less  chemical  dissociation  of  the  gases  due  to  a 
lower  static  temperature,  and  a less  severe  structures  requirement  as  a 
result  of  the  lower  static  pressures  and  temperatures . Heat  addition  in 
a supersonic  flow  may  lead  to  a strong  detonation  or  a weak  detonation 
wave  as  shown  in  Figure  (1)  which  represents  the  pressure-specific  volume 
relationships  between  initial  and  final  state.  The  Rayleigh  lines  repre- 
sent conditions  of  constant  mass  flow  with  varying  energy  exchange,  while 
the  Hugoniot  curves  represent  (for  a fixed  energy  input  and  initial  static 
temperature  and  pressure)  varying  mass  flow,  Mach  number,  total  tempera- 
tures, and  final  static  temperature  and  pressure.  The  weak  detonation 
mode  represents  a true  supersonic  combustion  process,  supersonic  flame 
speed,  which  does  not  involve  shock  waves  of  any  kind;  the  flows  of  both 
heated  and  unheated  gases  are  supersonic.  The  only  known  example  of  a 
thermal  initiated  weak  detonation  is  the  phenomena  of  “knock"  in  a piston 
engine  in  which  the  combustion  propagates  almost  at  an  infinite  speed 
throughout  the  gas  mixture.  In  the  strong  detonation,  a normal  shock 
precedes  the  combustion  process  so  that  energy  addition  takes  place  in 
the  flow  which  is  subsonic  relative  to  the  shock  wave.  The  Chapman- 
Jouguet  mode  is  the  strong  detonation  case  in  which  the  Mach  number  of 
the  flow  after  heat  addition  is  one  relative  to  the  shock  wave.  Whereas 
the  traveling  strong  detonation  waves  are  easily  obtained  and  have  been 
the  subject  of  many  studies,  stationary  strong  detonation  waves  resulting 
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from  chemical  energy  release  are  very  difficult  to  obtain  or  maintain 
because  of  difficulties  in  stabilizing  the  wave.  Standing  weak  detona- 
tions in  fuel-air  mixtures  are  difficult  to  achieve  because  of  the 
difficulty  in  obtaining  ignition.  Flame  speed  and  ignition  is  a strong 
function  of  the  static  gas  temperature  and  high  static  temperatures  are 
difficult  to  obtain  in  a supersonic  flow. 

The  final  pressure  and  volume  conditions  denoted  by  "x"  on  Figure 
(1)  can  be  reached  by  the  weak  detonation  as  discussed  above  or  by  a 
subsonic  diffusion  flame  propagating  in  a supersonic  stream.  The  weak 
and  strong  detonation  modes  can  occur  only  in  premixed  mixtures  of  fuel 
and  air.  Combustion  in  non-preniixed  gas  mixtures  depend  primarily  on 
the  rate  of  injection,  mixing,  diffusion  of  the  fuel,  and  chemical 
kinetics.  The  length  of  the  heat  addition  region  for  diffusion  flames 
will  be  much  longer  than  the  weak  detonation  because  of  the  processes  of 
mixing  and  diffusion  and  the  fact  that  the  normal  component  of  the  flame 
velocity  with  respect  to  the  incoming  mixture  does  not  propagate  at  a 
supersonic  speed  in  diffusion  flames  which  results  in  a small  flame 
angle  to  the  flow. 

There  is  a wealth  of  papers  dealing  with  the  problem  of  energy 
addition  to  a supersonic  flow  and  subsonic  diffusion  flames  propagating 
in  a supersonic  stream.  Only  a few  of  the  important  ones  will  be  refer- 
enced here. 

Weber  and  MacKay  [1]  showed  that  if  achieving  sonic  flow  at  the 
combustor  exit  was  not  possible  due  to  structural  or  fuel  limitations, 
then  maintaining  a shock  free  flow  during  energy  addition  is  desirable 
for  good  overall  engine  efficiency.  They  also  noted  that  the  thrust  and. 
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Figure  (1).  Hugoniot  and  Rayleigh  Curves  for  Heat  Addition. 
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in  some  cases,  the  engine  efficiency  can  be  increased  by  using  a 
combustor  with  an  increasing  flow  area.  There  exists  a combustor  tem- 
perature below  which  the  supersonic  combustion  ramjet  (SCRJ)  has  a better 
overall  engine  efficiency  compared  to  the  conventional  ramjet  (CRJ). 
Compared  with  a double  cone  inlet  CRJ,  the  isentropic  inlet  SCRJ  gives 
better  engine  efficiencies  above  a flight  Mach  number  of  5,  while  with  a 
wedge  inlet  the  SCRJ  is  more  efficient  above  M = 7.  They  noted  that 
friction  can  cause  significant  loss  in  total  pressure  if  the  combustor 
length  is  very  great;  thus  the  weak  detonation  wave,  which  would  have  a 
short  length  compared  to  the  subsonic  diffusion  flame,  would  be  advanta- 
geous here. 

Billig  [2]  describes  what  is  believed  to  be  the  first  demonstration 
of  stationary  supersonic  energy  addition  in  form  of  a subsonic  diffusion 
flame  propagating  in  a supersonic  stream.  He  noted  that  the  energy 
released  was  low  compared  to  the  maximum  theoretical  value  due  to  short 
residence  time  in  the  experimental  two-dimensional  combustor.  According 
to  this  investigation,  the  use  of  one-dimensional  equations  to  describe 
the  flow  is  possible  in  some  cases  but  cannot  be  adopted  "a  priori". 
Cookson,  Flanagan,  and  Penny  [3]  studied  subsonic  diffusion  flames 
propagating  in  a supersonic  stream  for  free  jets  and  enclosed  (duct) 
flow.  They  found  that  for  a free  jet  the  flame  was  easily  maintained, 
but  when  enclosed  in  a duct,  it  became  a strong  detonation  wave  unless 
the  duct  had  a certain  divergence  angle.  Too  large  a divergence  angle 
would  extinguish  the  flame.  They  used  a one-dimensional  analysis  for 
comparison  with  experiments  where  pressure  and  area  ratios  were  exponen- 
tial functions  of  distance  or  the  area  ratio  is  a function  of  the 
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pressure  ratio  (before  and  after  energy  addition)  raised  to  some  power. 
Their  pressure  measurements  were  within  5%  of  those  predicted  by  the 
theory;  the  theory  qualitatively  agreed  with  the  Mach  number  variation 
as  a function  of  total  temperature  variation. 

The  weak  detonation  wave  is  difficult  to  generate  in  the  laboratory 
with  combustible  mixtures  since  the  absence  of  heating  by  a shock  wave 
requires  heating  the  mixture  of  a supersonic  gas  flow  or  a photochemical 
process  to  obtain  ignition.  The  chemical  reactions  or  burning  must  take 
place  very  fast  (supersonic  flame  speed)  if  a weak  detonation  wave  is  to 
be  produced.  Heating  of  a supersonic  flow  of  a fuel-air  mixture  is 
virtually  impossible.  Pre-heating  a high  enthalpy  mixture  in  a stagna- 
tion reservoir  before  expansion  may  cause  premature  ignitions.  Thus,  an 
experimental  investigation  of  a thermally-initiated  weak  detonation  wave 
appears  to  be  very  difficult  in  the  laboratory.  Injection  of  the  fuel 
into  supersonic  flow  of  the  proper  total  temperature  and  Mach  number  will 
result  in  oblique  shocks  behind  which  subsonic  diffusion  flames  propagate 
releasing  a large  amount  of  the  energy  in  the  region  near  the  surface 
where  injection  takes  place. 

One  possible  approach  for  obtaining  a weak  detonation  is  by 
photochemical  ignition  of  a combustible  mixture.  Thermal  ignition  delay 
times  for  fuel-air  mixtures  are  rather  long  when  static  pressures  and 
temperatures  are  only  slightly  above  the  thermal  explosion  limit.  These 
long  times  can  result  in  very  long  combustor  lengths,  thus  if  photo- 
chemical ignition  is  possible,  then  the  ignition  times  may  be  substan- 
tially less  resulting  in  shorter  and  lighter  weight  combustors. 
Experiments  by  Levy,  Cerkanowicz,  and  McAlevy,  III  [4]  demonstrated 
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that  ical  ignition  delay  times  are  possible;  they 

investigated  flash  photolysis  of  hydrogen-air  and  methane-air  mixtures 
in  the  vacuum  ultraviole"  from  1350  to  1700  Angstroms.  Rice  [5]  inves- 
tigated the  formation  of  weak  detonation  wave  by  exposing  hydrogen- 
■ me  and  chlorine  sensitized  oxygen  hydrogen  mixtures  to  visible 
and  near  ultraviolet  light.  A shock  tube  was  used  to  generate  the  gas 
flow  and  a Xenon  flash  lamp  provided  the  radiation.  Experimental  mea- 
surements and  theoretical  calculations  showed  that  the  chemical  reactions 
resulting  from  photochemical  initiation  are  fast  enough  to  form  a weak 
detonation.  However,  to  sustain  the  combustion  of  the  system  studied, 
very  high  power  levels  on  the  order  of  hundreds  of  kilowatts  are  nec- 
essary for  the  photon  source. 

The  purpose  of  this  dissertation  was  to  generate  in  constant  and 
varying  area  a weak  detonation  wave  by  radio  frequency  (R.F.)  heating  of 
a supersonic  inert  gas  and  by  mass  addition  to  a supersonic  water  stream. 
To  establish  the  supers'  ve,  an  R.F.  coil  is  wrapped  around 

an  electrically  non-conducting  duct  through  which  gas  flows.  The  coil  is 
connected  to  an  R.F.  power  supply.  The  length  of  the  coil  will  represent 
the  wave  region  since  the  heatii  is  accomplished  by  electrons  and  con- 
sequently will  be  quite  rapid.  Thus  one  could  measure  the  flow  properties 
of  a supersonic  flow  through  the  duct  and  compare  them  to  the  equations 
for  a weak  detonation.  Since  heating  by  an  R.F.  coil  restricts  the  major 
energy  input  to  a region  about  one-half  the  distance  from  the  wall  to  the 
center  of  the  duct  because  of  space  charge  effects,  a variation  of  tem- 
perature with  duct  radius  can  be  studied  so  that  a good  check  can  be 
made  on  the  ideal  qua:,i-one-dimensional  equations  which  are  based  on 
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average  properties  across  the  duct  radius.  R.F.  heating  of  a gas  is 
relatively  inefficient,  only  about  10%  of  the  power  generated  can  be 
transferred  to  the  gas;  hence,  for  an  actual  engine  this  would  have  to 
be  greatly  increased  before  R.F.  heating  can  be  considered  as  a practical 
form  of  energy  in  powerplants. 

Oppenheim  and  Teamen  extended  the  hydraulic  analogy  theory  to 
include  the  energy  addition  case  for  comparison  of  water  and  gas  flows. 

In  the  basic  analogy  water  behaves  like  a gas  with  a specific  heat  ratio 
(y)  equal  to  two  and  gas  pressure  ratios  are  represented  by  water  height 
ratios  squared.  The  energy  addition  analogy  is  qualitative  and  not 
quantitative.  Oppenheim  [6]  and  Teamen  [7]  experimentally  investigated 
the  case  of  the  strong  detonation  wave  (hydraulic  jump  created  by  mass 
addition  to  the  main  stream  flow)  but  did  not  experimentally  check  the 
weak  detonation  mode.  In  this  dissertation  the  simulation  of  a weak 
detonation  wave  by  mass  addition  to  a supersonic  water  flow  was  attempted 
to  determine  if  it  would  qualitatively  predict  the  phenomena  obtained  in 
the  gas  flow  tests.  This  effort  would  validate  the  use  of  the  hydraulic 
analogy  to  obtain  an  insight  into  energy  addition  to  a flowing  gas  for 
any  velocity. 
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SECTION  III 
THEORETICAL  ANALYSIS 


Flow  Equations 

The  quasi-one-dimensional  theory  for  gas  flow  is  based  on  average 
properties  (with  respect  to  radius,  =J  MRdR/ J RdR)  across  the 
radius  of  an  axisymmetric  duct.  From  Shapiro  [8],  one  may  find  the  spe- 
cial conditions  for  energy  addition  based  on  constant  area,  constant 
Mach  number,  constant  pressure,  etc.  However,  for  the  general  case  of 
energy  addition  simultaneous  with  area  variation,  friction,  and  mass 
exchange,  a different  approach  must  be  taken.  For  a comparison  with 
experimental  results  which  are  based  on  energy  exchange  in  varying  area 
ducts,  it  is  necessary  to  have  a theory  which  includes  these  effects. 

Two  approaches  were  considered.  A rigorous  analysis  retaining  all 
terms  of  the  fundamental  equations  is  very  tedious  and  entails  great 
difficulties  in  working  out  a computer  program.  This  first  approach 
starts  with  the  Navier  Stokes  equations  and  makes  no  simplifying  assump- 
tions. To  solve  the  equations,  initial  profiles  (along  with  radial 
derivatives)  across  the  radius  for  pressure,  temperature,  and  velocity 
are  assumed,  and  all  axial  derivatives  are  written  in  finite  difference 
notation.  The  energy  input  term  is  given  as  a function  of  radius  and 
axial  position.  The  finite  difference  terms  are  written  based  on  one 
point  in  front  of  present  location  and  up  to  four  points  behind  it. 

This  method  allows  one  to  solve  for  the  AU,  AT,  etc.  that  are  required 
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for  going  from  the  present  location  to  the  next  one.  The  solution  is 
obtained  at  each  point  across  the  radius  by  varying  the  pressure  incre- 
ment until  an  agreement  is  reached  between  the  radial  velocity  component 
in  the  continuity  equation  and  the  radial  velocity  component  in  the 
radial  momentum  equation.  Once  this  calculation  is  accomplished  across 
the  radius,  the  values  of  the  properties  are  known  and  their  first  and 
second  derivatives  can  be  calculated  by  a SPLINE  computer  routine  which 
fits  a curve  through  a set  of  given  data  points.  After  the  new  deriva- 
tives are  known,  the  process  is  repeated  for  the  next  axial  increment. 
The  difficulty  of  the  above  analyses  lies  in  providing  a realistic 
initial  profile  (particularly  near  the  wall)  for  the  axial  and  radial 
velocities.  Too  large  an  increment  in  velocity  was  usually  calculated 
near  the  wall  or  a solution  was  not  obtainable.  When  a solution  was 
obtained  across  the  radius,  calculation  of  the  second  derivative  of 
velocity  with  respect  to  radius  resulted  in  values  with  alternating 
signs  thus  preventing  continuation  of  the  solution.  Because  of  these 
difficulties,  this  effort  has  been  discontinued  until  a better  method 
for  calculating  second  derivatives  can  be  developed. 

In  a second  approach  the  problem  has  been  analyzed  successfully  by 
a procedure  which  is  an  extension  of  Shapiro's  [8]  method.  Following 
Shapiro,  the  Navier  Stokes  equations  are  used  for  the  case  that  the 
coefficient  of  viscosity  is  constant,  contributions  to  the  shear  stress 
by  ?u/5x  is  negligibly  small,  the  flow  is  axisymmetric,  average 
properties  (P)  defined  by  (P)  = JTR(P)RdR/jfR  RdR  are  used,  and  a fric- 
tion factor  f is  defined.  For  this  case  it  can  be  shown  that: 
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Equation  of  State  (4) 


D is  the  diameter  of  the  duct  and  y is  the  injected  mass  flow  velocity 
component  parallel  to  the  main  stream  velocity  divided  by  the  main 
stream  velocity.  Defining  a conical  flow  area  we  can  write: 


P = m/p  1-  D, 


el  A _ 2 dP  _ 2m  Jv 
A 0 rrw  + D,  (5) 


where  m = 2tan9  and  0 is  the  divergence  angle  of  the  duct.  In  the 
following  work,  the  average  sign  (-)  will  be  dropped  from  the  equations. 

Multiplying  the  momentum  equation  (2)  by  dx  and  dividing  by  pressure 
p and  using  equation  (1)  we  get: 
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dul  = dM%  Jl°  _ M 

u?  Ma  T°  If  ~ ^ M* 


(8) 


To  extend  Shapiro's  analysis,  the  following  two  assumptions  are 


made: 


These  two  assumptions  will  permit  the  analysis  of  a flow  in  which  total 
temperature  change,  area  change,  and  mass  flow  exchange  can  occur 
simul taneously . 

Substituting  equations  (1),  (3),  (4),  (8),  (9),  and  (10)  into  (6) 
and  rearranging  results  in: 
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To  determine  the  proper  integral  to  use  for  evaluating  the  integral  of 


- dM 


M2(a  + t>M2+cM‘') 


on  the  left  side  of  equation  (11),  it  is  nec- 


essary  to  examine  the  quantity  —b  , or: 


__  [y lfA  + N+  anf+-  R*-C2.~N-2.n).'tl .( A.  +N+zn) 
. Vm  J 2 \srn  / 
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Taking  the  first  derivative  of  q with  respect  to  <j>/2m,  and  setting  equal 
to  zero  gives: 

ITT  = - *r  ] (,3) 

Also  setting  q = 0 and  solving  for  $/2m,  we  get: 

m = -[T7-(N+2")+ 

Taking  the  second  derivative  of  q with  respect  to  <f>/2m,  we  get: 


i &lm)1 


- ~zr 


Hence,  since  the  function  is  concave  downward  and  the  maximum  occurs  at 
the  same  <|>/2m  value  where  q = 0,  q is  always  £ zero.  This  condition 
determines  the  type  of  integral  to  use  from  the  integral  tables. 

Now  integrating  both  sides  of  equation  (11)  and  combining  terms, 
we  final ly  arrive  at: 
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where 


y = 


(X-f-l)N  +2[yC %m-H)  - 0-n)] 


Equation  (16)  allows  us  to  calculate  the  other  properties  for  heat 
addition  with  simultaneous  varying  area,  mass  flow  exchange,  and  fric- 
tion by  utilizing  the  equations: 
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which  result  in: 
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(21) 
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With  these  equations  [(16)  through  (23)],  using  n or  N as  variable 
parameters,  and  including  friction  by  means  of  the  friction  factor  f, 
one  can  analyze  flows  in  axisymmetric  ducts  with  area  change,  energy 
exchange,  and  variable  mass  flow.  The  assumptions  (9)  and  (10)  are 
reasonable  since  eight  special  cases  satisfy  this  relation: 

1.  Zero  energy  exchange  N = n = 0.0 

2.  Constant  area  duct  N -*•  °°  n -*  °» 

3.  Constant  Mach  number 

N - id-A-yfw?  +2(i-  V)n]  M? 

vm;  + 1 
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f » friction  factor 
6 - di vergence  angle 
y - injection  velocity/flow  velocity 

4.  Constant  total  pressure 


N 


f 


5.  Constant  static  pressure  (M  0.5) 


N 


I -n  - 


6.  Constant  density  and  Mach  number 


w = - r 2 a-  yj n ) 

t / ) 


7.  Constant  temperature  and  Mach  number 

K - o.o 


-f 


- 0-n) 


- en 0-  *) 


8.  Constant  pressure  and  Mach  number 

n = —fit -2  0-  >)n 

w? r+2°-J,)'n 
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The  special  case  (6)  (constant  density  and  Mach  number)  will  be 
demonstrated  here.  Starting  with  equation  (6)  we  can  write: 


(24) 
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gr  + io-wllr 


JL”=  N J± 
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^ J 2 L 2m 

Therefore,  for  constant  density  and  Mach  number: 

N = - ¥(l£  + *c'-^n)/(f1^i)  ,2£ 

It  should  be  mentioned  that  the  constant  total  pressure  condition 
(4)  cannot  be  achieved  if  the  energy  exchange  is  dependent  on  friction 
by  Reynolds  analogy;  and,  condition  (4)  requires  energy  removal  from 
the  flow  to  counteract  the  effect  of  friction.  Condition  (4)  does  not 
hold  for  heat  addition.  Total  temperature  based  on  Reynolds  analogy 


can  be  written: 


5 = -i) 


wall  temperature  (29) 


From  special  condition  (4)  for  zero  mass  transfer,  we  can  write  the 


following: 
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Equation  (30)  can  also  be  directly  obtained  from  the  momentum  equation 
(2)  by  the  appropriate  substitutions  resulting  in: 


ip  y/vf 

P*  2 
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r°  2 

which  for  constant  total  pressure  results  in  equation  (30).  Equation 
(30)  is  different  from  equation  (29),  thus  the  cooling  of  the  flow  must 
be  achieved  by  an  ideal  process  which  tends  to  lower  the  gas  temperature 
and  not  by  heat  transfer  across  the  duct  wall. 

If  LN  (n)  is  zero,  there  is  no  friction,  and  N -*  ®,  then  equation 
(16)  reduces  to: 
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Now,  if  this  is  correct,  then  for  a normal  shock  (q/CpT^  = 0.0)  the 

right-hand  side  of  equation  (31)  reduces  to  1.0.  To  check  this  result. 
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let  =4,  the  Mach  number  behind  the  shock  = 0.435,  and  y = 1.4. 
Substituting  into  equation  ( 31 ) gives: 
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(33) 


Hence,  we  have  demonstrated  that  equation  (16)  is  correct  and  equation 
(31)  holds  for  constant  area  and  no  friction. 

There  are  probably  some  values  of  N and  n and  for  which  the 
conditions  of  constant  pressure,  constant  temperature,  or  constant 
density  are  approximated;  however,  it  will  be  more  convenient  to  start 
with  the  basic  equations  to  evaluate  these  conditions  retaining  only  the 
expression  m^/rfi-j  =[A2^Aijn‘  Below'  the  necessary  equations  for  calculat- 
ing the  constant  conditions  mentioned  above  are  given. 

Constant  Pressure 


2 

I - n + P/*irr\ 


(34) 


Ji 

r, 


(2(1- n-  t/nm)  \ 
M.  1 V / 


Constant  Temperature 


(35) 


(36) 
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C n-/) 


For  friction 


A - fJL\ 

a,  l*w 


Constant  Density 

— 

A - jMi_ 
A.  ~ 


*2  (/ 


( i-n ) 


'-(n t ife 


- Hm(i~n)\ 
24>0-n)  / 


~£  0 1/m(n-i))M,?f  I 

f 0 + + 1 


Hrr\  tn-i) 


-H + 1 


Now  a few  special  conditions  that  can  occur  in  equation  (16)  will 
be  examined.  Consider  the  exponent 


T = 


0*0/\/  + Z[y(i-+  <p/dm)-u  n~] 


and  the  expression 


i + r' 


i * *r  M? 


20-n)-/v/-yr4>/Pmt^)Alf 


= 7 


When  the  denominator  of  y is  0,  then  ^ ->  « and 

N = -s  »('  + &) -O-nilyci/i-o 


y 7 I 
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A , 

A. 


finite  when  n = 1 


(45) 


ft. 


0 when  n < 1 


Writing  • 1 and  solving  for  4>/2m , it  is  found  that: 


<t> 


> 

2m  ^ 


2*  * 

also,  from  (43)  we  can  write: 

> y-»-/ 


*+l  N+  -^r1-  - (i-  n) 


■N  + 


i-n 


- i 


(46) 


(47) 


(48) 


2m  ^ 2*  Y 

The  last  two  equations  are  equal  for  n = 0;  hence,  when  ^ the 
equation  then  reduces  to  n'  = (1)  =1  and  a solution  exists  to  equation 

(16).  Next,  if  the  terms  in  equation  (16)  of  cZ(l-h)-N~ll/l/\2(,2m  . 

a.  + bM2  + c Mg  and  + b t | (-fy  are  set  equal  to  zero,  it  is 

found  that: 

2 _ 2Q-n)-A/ 

^2  y (4/zm  + N ) (49) 

for  each  term  and  likewise  for  the  corresponding  terms: 

l _ 2 Q-n)~  /V 

y (<l>/2/ntA') 

For  the  constant  Mach  number  case  M2  = M^ ; hence,  solving  for  N one 
arrives  at:  ^ 

N - zO-n)  - y + 2Q-*)n) 

*/n2  + / 


(50) 


(51) 
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which  is  the  third  special  case  listed  earlier.  If  N = 0 = n,  equation 

2 

(51)  gives  tane  = yfM  /2,  which  is  the  case  of  constant  Mach  number  for 
friction  as  presented  in  Shapiro,  Vol . I [18].  From  the  above  it  is 
seen  that  for  a solution  to  exist  for  constant  M then: 


2 6-n)  - N >0.o 


(52) 


In  equation  (16)  another  term  that  needs  to  be  examined  is  the 
exponent  1 /2a  = 1/(2  [2  (1-n)  - N])  when  the  denominator  equals  zero; 
2 (1-n)  - N = 0 or  N = 2 (1-n).  When  a is  zero,  b/Y^T  = 1 or 
b = -y  ($/2m  + N),  c = (y-1 )/2  b,  and  equation  (16)  reduces  to 


Aj 

A. 


_M? 

M2 


i * 

W2  + 2.  m 2 


tlMf  + 


+ I 

' 2.  • 2.  2 

*-/  alY 

/ -f-  Z ^2 


L/  + 


Mj 


or 


A 

A 


= 0) 


It  TM2 


JVC 


1 !/- 


<?.  o 


(53) 


(54) 


Az  _ 

'it  ‘ 

Mi 

A,  ' 

1 J 

• 

_ m 

2. 


(55) 

(for  a = 2 (1-n)  - N = 0.0) 


Hence,  a solution  exists  for  this  condition. 

Equation  (16)  has  been  discussed  and  analyzed  in  some  detail  so  it 
is  obvious  that  a large  range  of  flow  conditions  can  be  investigated 
with  it.  Equation  (16)  is  used  for  comparing  theory  with  the 
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experimental  results  in  Tables  2-10.  Some  examples  of  flow  utilizing 
equation  (16)  are  shown  in  Figures  (2)-(9).  These  figures  represent 
flow  with  area  change,  mass  change,  and  total  temperature  change 
occurring  simultaneously.  For  LN  (or  n)  equal  to  1000,  this  represents 
the  constant  area  case.  Figures  (10)-(14)  represent  area  change 
(decreasing)  simultaneous  with  total  temperature  change  with  no  fric- 
tion. For  N equal  1000,  the  area  ratio  equals  one.  The  fact  that 
equation  (16)  reduces  to  the  constant  area  for  LN  (n)  or  N + » 
can  be  checked  by  choosing  values  off  the  above  figures  and  comparing 
them  against  the  constant  area  case  in  Shapiro  [8].  The  figures  were 
obtained  by  a computer  program  which  was  written  to  analyze  equations 
(16)  through  (21). 

As  an  example  of  how  to  use  the  curves,  assume  you  have  a duct 

with  an  area  ratio  of  Ag/A-j  = 2.0,  a mass  addition  ratio  of  1.6,  and 

-0  5 

the  total  temperature  ratio  proportional  to  (Ag/A-|)  ' ’ while  M-j  = 4.0. 

From  Figure  (2),  this  corresponds  to  n or  LN  = 0.7,  which  is  one  of  the 
parametric  curves  based  on  various  n values;  m^/m^  = (A^/A^)11.  From 
Figures  (3)  and  (4)  it  is  found  that  Mg  = 1.5.  Figures  (5)-(9)  then 
give  the  following  values:  Tg/T°  = 0.7,  Tg/T-j  = 2.1,  P^/P-j  = 3.75, 

Pg/P°  = 0.075,  and  Vg/V-j  = 0.54.  If  one  is  interested  in  an  initial 
Mach  number  other  than  M^  = 4.0,  then  it  is  necessary  to  take  the  ratio 
of  the  y axis  values  for  the  desired  initial  and  final  Mach  numbers. 
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NOTATION  FOR  FIGURES  (2)  THROUGH  (14) 

G -*■  y 

Ml  * M1 , PI  - P] , etc. 

LN  - n 

Ml  ■*  m-j  mass  flow  rate 

PH  -*■  friction  parameter  (f/tana) 

N = O.OC  ■*  constant  total  temperature 
N = 1000  -*■  constant  area 
n = LN  = 1000  -*■  constant  area 
n = LN  = 0.0  -*•  0.0  mass  flow 

Mass  flow  variation  assumed  to  occur  by  parallel  injection  or  removal 
of  mass. 
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Water  Table  Theory 

It  was  noted  in  the  introduction  that  there  exists  an  analogy 
between  gas  flow  and  water  flow,  both  subsonic  and  supersonic,  which 
makes  it  possible  to  study  aerodynamic  problems  on  a water  table.  The 
analogy  will  now  be  presented  in  a mathematical  form.  See  Ref.  [6]  and 
Ref.  [7],  which  discuss  the  strong  detonation  wave  analogy. 

For  steady-state  water  flow  in  a constant  area  channel  we  can 
write:  (h  = water  depth;  w = channel  width)  (mass  injection  normal  to 
the  flow) 

• i2  • 2 

mv+g;fw_h  = m(v,  i-  qf  WJt,  (56) 

2 Z 


mL 

sw 


(57) 


Vh 


+ 5 


H 


(58) 


(59) 
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Where  h = = / + 

qs  = injected  mass 
For  gas  flow  we  can  write: 

P P S Vz  = P,  t V? 


% 


l 


(60) 


2 2 

By  comparing  (59)  and  (60),  we  see  an  analogy:  gh  /2  p,  X /h  v 

(specific  volume).  Pressure  is  proportional  to  the  water  depth  squared. 

2 2 
Let  gh  /2  = y be  equivalent  to  pressure.  Let  x /h  = x be  equivalent  to 

specific  volume.  Thus,  equation  (59)  becomes: 


"»!  *v  = *.  »f  + 

or  «» = r t- f_  r x 1 

/ - y/y,  ' 

* l a J [ *.  J 

.*/*,  - / _ 

(61) 


which  is  the  Rayleigh  Line  analogy.  The  slope  of  the  Rayleigh  Line  is 
tana  = dy/dx  = q^  = -h^  , or  ^tana  = h^V-j,  which  is  similar  to  the 


gas  case  where  ytana  = y M] . 

Defining  F as  the  ratio  of  velocity  to  gravity  wave  speed,  we  have: 


r _ VL  _ ±/iL  - I I 

F - c ' VsiT  “ * ! 


ay 


(62) 


and  for  a gas 

m = i 


m I 

- 


mr  = * i »w 

Thus,  by  analogy  m/A  + q/x  and  y ■>  2. 


FT 


(63) 
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I 


It  is  apparent  that  water  flow  is  analogous  to  a gas  flow  having  a 
ratio  of  specific  heats  equal  to  2.  The  Rayleigh  analog  can  also  be 
written: 


i - v/y. 


x/x,  - I 

and  the  gas  case  as 

) ~ P/f? 


1_  = 

V*,  ' 


= vm; 


}/y  - < 

From  the  definition  of  entropy 


<JS  + R jhL 

J n V 


AS 


= ^ t,  + 


(*-i) 


A 

p 


AS/Cv 


we  can  write  for  a gas  that: 


y A S/Cy, 

pv>  = ffv/c 


Thus  to  get  an  analogy  for  water,  we  can  calculate: 

y x2  = Z3  vy 


2 - S' 


y*  = 


/ 


A 

S; 


-.V 


_ 

£ 


I + % 


x,  - X 


v - y, 


•i 


(64) 


(65) 


(66) 


(67) 


(68) 


For  the  initial  point  y^2  - g/2 ; hence,  we  can  write  the  Hugoniot 
analog  as: 
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l ~ */X, 

y/y.  - i 


¥ 


kt  -i  2- 

(>t)  ’ C y/y,  - / ) 

TVy,  • (yA,  -i)  +2  e* 


(69) 


where 


/ f /- 


(%,  -<-i e’XWr,  ■ C Vy, -i )*!(*) 
Wx  ( vy,  - if 


. -1^ 


(70) 


This  Hugoniot  curve  for  mass  addition  differs  from  the  gas  heat 
addition  case  since  it  passes  through  the  point  y/y-j  = x/x-j  = 1.0  and 
the  gas  case  doesn't.  The  above  analysis  shows  that  an  analog  exists 
between  mass  addition  to  a water  flow  and  heat  addition  to  a gas  flow. 
Although  one  cannot  get  a quantitative  analysis  for  the  gas  flow  from 
the  water  table,  it  is  possible  to  gain  a qualitative  insight  of  what  to 
expect  in  a similar  gas  flow  situation.  The  nozzles  fabricated  for  the 
water  table  were  based  on  a sharp-edged  throat  and  a gamma  of  2.0 
utilizing  the  two-dimensional  approach  outlined  in  Shapiro,  Vol . I [8]. 


Electric  Probe  Theory 

To  analyze  the  results  of  the  R.F.  energy  addition  experiments,  it 
is  necessary  to  make  total  temperature  measurements  across  the  duct  exit 
with  a total  temperature  probe.  However,  since  there  is  ionization  of 
the  gas,  it  is  necessary  to  account  for  the  energy  released  in  the  probe 
by  electron  recombination.  An  electric  double  probe  was  used  to  deter- 
mine the  ionization  level  so  that  a correction  could  be  made  to  the 
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measured  total  temperature  in  order  to  account  for  electron-ion 
recombination  and  the  electron  temperature.  The  probe  theory  will  now 
be  discussed  in  some  detail. 

The  original  single  probe  theory  as  developed  by  Langmuir  [9] 
applied  to  plasmas  (ionized  gases)  where  the  mean  free  path  of  the  par- 
ticles was  much  larger  than  the  radius  of  the  probe.  According  to  this 
technique,  a voltage  is  applied  from  the  probe  to  one  of  the  discharge 
electrodes  and  the  resulting  current  is  measured  as  the  voltage  sweeps 
from  negative  to  positive.  Langmuir  showed  how  the  electron  temperature 
can  be  determined  from  the  slope  of  the  electron  current  vs.  voltage 
curve  and  that  the  electron  density  is  derived  from  the  current  magni- 
tude and  electron  temperature.  See  Figure  (15)  for  a sketch  of  the  probe 
characteristic.  The  negative  voltage  part  of  the  curve  is  the  ion  collec- 
tion part  of  the  characteristic;  this  is  the  part  of  the  curve  that  is 
relevant  to  the  double  probe  theory. 

Johnson  and  Mai  ter  [10]  developed  a theory  for  the  use  of  a double 
probe  for  measuring  electron  temperature.  They  pointed  out  that  the 
single  probe  method  (SPM)  is  valid  only  if  the  cathode  plasma  potential, 
electron  temperature,  and  random  electron  current  do  not  change  with  the 
cathode  probe  potential.  However,  as  the  electron  current  (i  ) increases, 
one  often  observes  changes  in  the  discharge  patterns,  especially  if  i 
is  an  appreciable  fraction  of  the  main  discharge  current.  The  probe  is 
then  disturbing  the  very  quantities  it  is  intended  to  determine.  As  can 
be  shown  for  the  double  probe  method  (DPM),  the  total  current  to  the  probe 
can  never  exceed  the  positive  ion  current  to  them;  and,  since  this  cur- 
rent is  hundreds  of  times  smaller  than  the  electron  current  for  the  SPM, 
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Vp  - Probe  Potential  w.r.t.  Arbitrary  Reference  Point 
V^.  - Floating  Potential,  Electron  and  Ion  Currents  Are  Equal 

Vs  - Space  Potential,  Probe  at  Same  Potential  as  Plasma 
A-B  - Electron  Collection 

B-C  - Ion  Collection,  Part  Relevant  to  the  Double  Probe  Analysis 


Figure  15.  Schematic  of  the  Langmuir  Probe  Characteristic  Curve. 


»6 
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the  DPM  appears  to  have  an  advantage  in  this  respect.  Also,  in  an 
electrodeless  type  discharge  as  used  in  the  experiments  to  be  discussed 
later  where  there  is  no  reference  electrode  for  the  SPM,  it  is  advanta- 
geous to  use  the  DPM. 

The  equation  developed  by  Johnson  and  Mai  ter  for  the  electron 
temperature  is: 


where 


e = electronic  charge 

k = Boltzman  constant 


d Vd /d i d = the  equivalent  resistance  or  slope  of  the  current  vs. 
voltage  curve  at  the  origin. 

See  Figures  (16)  and  (17)  for  diagrams  illustrating  the  above  quantities. 


For  the  case  where  the  probes  are  of  equal  area,  random  current  density 
at  both  probes  are  equal,  the  constant  emf  between  the  probes  is  zero, 
and  the  space  potential  near  both  probes  is  the  same,  then  equation  (71) 
reduces  to 


If  the  current  magnitudes  are  not  the  same  but  the  curve  passes  near  or 


through  the  zero  point,  then  a good  approximation  can  be  made  by  writing: 
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Ideal  Case 


Figure  16.  Voltage-Current  Character! Stic  of  the  Double  Probe. 


Figure  17.  General  Potential  Diagram  for  the  DPM. 


jQ1  = Electron  Space  Current  Near  Probe  #1 
jo  = Electron  Space  Current  Near  Probe  £2 
= Probe  to  Plasma  Potential  Probe  #1 
V£  = Probe  to  Plasma  Potential  Probe  #2 

= Contact  EMF  Between  Probes  (=  0 for  Symmetrical  Probes) 
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r,  = 


_ jl  L* ' Lp* 


L,  ¥ L, p 
• i 


J Lj 


(73) 


= o 


For  the  nearly  symmetrical  probe  characteristic,  it  can  be  shown  that: 


r el/, 
L 


(74) 


The  traces  made  in  the  experiments  closely  followed  this  relationship, 
see  Figure  (18)  for  a typical  double  probe  trace  from  the  R.F.  experi- 
ments. Johnson  and  Malter  noted  in  their  article  that  only  the  electrons 
sampled  in  the  DPM  are  the  top  few  percent  of  the  distribution  and  thus 
the  electron  temperature  determined  will  be  somewhat  higher  than  the 
average  Maxwellian  temperature.  The  electron  mean  free  path  in  the 
experiments  was  about  ten  times  the  probe  diameter  and  Langmuir's  theory 
for  electron  collection  would  be  approximated. 

Mow  that  a method  is  available  for  determining  the  electron 
temperature,  it  is  necessary  to  find  an  equation  which  will  give  the  ion 
density  utilizing  Tg.  For  the  pressures  involved  in  the  experiments, 
it  was  not  possible  to  use  the  collision-free  Langmuir  model  for  ion 
collection.  For  density  measurements  at  high  pressures,  the  papers  of 
the  following  authors  were  consulted:  Waymouth  [11],  Cohen  [12],  and 

Su  and  Lam  [13].  These  investigators  analyzed  spherical  probes  that 
collected  current  in  dense  plasmas  where  the  mean  free  path  is  much 
smaller  than  or  on  the  order  of  the  probe  radius.  If  one  takes  their 
final  equations  and  reduces  them  for  ion  collection  at  high  negative 
potentials  and  Te  >>  T;,  they  all  approach  the  same  form  as  will  be 
shown. 
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Waymouth  based  his  derivation  on  the  assumptions  that  Xc/rp  < 10, 
rp/100  Ap  > 1.0,  and  *c/Ap  > 10,  where  Ac  is  the  particle  mean  free  path, 
Ap  is  the  Debye  length,  and  rp  is  the  probe  radius.  The  potential  of  any 
chosen  ion  or  electron  in  a plasma  drops  off  as  q/r-e"r,/  ^ instead  of  q/r 
as  in  free  space.  The  field  of  the  chosen  ion  can  be  regarded  as  being 


effectively  screened  off  at  distances  beyond  Ap  where 


Su  and  Lam,  along  with  Cohen,  assumed  that  the  probe  radius  is  large 
compared  with  the  Debye  length.  Su  and  Lam  assumed  very  large  probe 
potentials  and  very  large  T /T.j , while  Cohen  kept  them  finite  in  his 
analysis;  for  Cohen  Ac  <_  10  Ap.  Figure  (19)  shows  the  regions  of  applica- 
bility based  on  the  assumptions.  On  the  figure  are  points  calculated 
from  the  experiment.  Since  all  theories  reduced  to  essentially  the  same 
form  for  T /T.  >>  1 and  large  negative  potentials,  any  data  points  fall- 
ing within  the  Waymouth  and  Cohen  regions  can  be  analyzed  by  equation 
(97).  Those  few  points  falling  outside  the  region  were  near  the  center 
of  the  duct  where  the  electron  density  was  low  and  the  total  temperature 
increase  was  quite  small;  hence,  any  errors  here  would  have  a small 
effect  on  the  average  total  temperature  calculation.  See  Table  11  for 
indications  of  the  magnitude  of  ion  number  density  and  electron 
temperature. 

First  consider  Su  and  Lam: 


Ti  - ion  temperature  (76) 
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Figure  19.  Domain  of  Validity  for  Probe  Theories  with  Data  Points  Shown 
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<£ 

e+l 


C 


c - constant  order  (one)  (77) 


l*  = *ir I;  n.  m;  kT; 


• - ion  mobility 
nQ  - ion  density 


(78) 


I; 


C 
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Finally  from  Cohen 
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where  pp  is  the  radius  of  the  probe  divided  by  the  Debye  length  and 
X = e(J+-J_ )/(£ (J+)+J_<)  • X 1 as  the  potential  becomes  very  negative 
and  causes  J to  go  to  zero.  J+  can  be  written: 
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where  D is  the  diffusion  coefficient,  and 
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where  ui  is  the  ion  mobility.  Utilizing  the  tables  in  Cohen,  equations 
(88)  and  (89)  can  be  written: 
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J*  =K-)('+i Tf)s'»»('^)  (92) 

= ('+  6 ) ('  + 5T?  ) = '^(/f  -g-  ) (93) 

or 

= c ( / t 7e  /7t-  ) (94) 

where  C is  the  order  one;  thus,  we  can  write: 

= c ■ c V7T  rp  mi  Oo  (,  4.  re/rt-  ) (95) 

and  for  T >>  T . : 
e i 

I+  = C ^ TT  Tp  M-  K Te  (96) 

This  is  the  same  form  as  equations  (81)  and  (88). 

Thus  the  equation  used  for  determining  the  ionization  level  was 
derived  from: 
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where  Tg  is  determined  as  shown  earlier. 
The  ion  mobility  u.  is  determined  by 
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where  piQ  refers  to  the  density  at  standard  conditions  and  y..  is  the 
mobility  for  standard  conditions  which  is  a known  constant. 
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Mi  = Ml 
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Hence,  we  can  write: 


T = H>7rc, 

+•  ^ h 
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and  since  nQ/n  = a,  the  degree  of  ionization,  we  have: 


<*  = 2 G>1  X 10 


wrp  H,  .o  re 


I+  •+  saturation  ion  current  (101) 
measured  from  I vs.  V 


This  is  the  equation  used  to  determine  a for  correcting  the  total 
temperature  measurement.  The  radius  of  the  probe  r is  for  a spherical 
probe;  the  probes  used  in  the  experiments  were  cylindrical.  Waymouth 
[11]  suggested,  in  light  of  the  physical  mechanism  he  considered,  that 
two  probes  that  draw  the  same  current  from  a given  plasma  will  cause  the 
same  perturbation;  and  since  probes  draw  current  in  proportion  to  their 
areas,  one  could  say  the  cylindrical  probe  has  an  effective  radius  equal 
to  the  radius  of  the  sphere  with  the  same  surface  area: 


rP  = C±  n,  ) 

c-f  f 


for  wire  probes, 


where  r^  is  cylindrical  radius  and  t,  is  length  of  the  wire. 

r is  the  radius  used  in  equation  [101]. 
peff 

Next  we  will  derive  the  equation  for  correcting  the  measured  total 
temperature  for  ionization. 


•jr  1 + h + ev  = h 
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EV  ->  Energy  for  ionization  of  argon  atom. 

Thus  assuming  complete  recombination,  one  may  write: 
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This  is  the  equation  used  for  correcting  the  total  temperature 

4 

measurement,  where  EV/Cptn^  for  Argon  is  equal  to  6.9  x 10  °K. 

With  regard  to  the  errors  arising  from  the  electrostatic  probe 
measurements,  the  comments  of  two  authors  are  cited.  Danikov  and 
Gruglyakov  [14]  indicate  that  at  pressures  on  the  order  of  1 mm  Hg,  the 
electron  density  determined  by  electrostatic  probe  measurements  could  be 
low  by  a factor  of  two.  Schulz  and  Brown  [15]  indicate  that  plasma  num- 
ber densities,  using  single  or  double  electrostatic  probes,  can  be 
determined  within  a factor  of  two  up  to  pressures  of  6 mm  Hg.  Thus, 
based  on  the  above  observations,  the  error  in  measured  electron  number 
density  or  a in  equation  (108)  should  be  on  the  order  of  two.  However, 
the  constants  dropped  from  the  Su  & Lam  and  Cohen  results  were  on  the 
order  of  1.1  and  1.2,  which  would  indicate  errors  up  to  20%  if  one 
ignored  the  higher  than  normal  electron  temperature  mentioned  earlier  in 
Johnson  and  Walter's  work. 

The  assumption  that  the  recombination  is  100  percent  in  the  total 
probe  and  all  the  energy  is  converted  into  thermal  energy  of  the  parti- 
cles is  a reasonable  one  based  on  the  results  of  a paper  by  Chen  [18]. 
From  his  work,  it  is  found  that  during  recombination  for  Argon  at  the 
electron  densities  found  in  the  experiments,  about  90  percent  of  energy 
released  is  transferred  to  the  electron  gas  and  10  percent  is  in  the  form 
of  radiation.  Tnus,  inside  the  total  temperature  probe,  most  of  the 
radiation  can  be  trapped  by  the  surrounding  walls  and  eventually  trans- 
ferred to  the  gas  molecules  through  collisions.  By  collisions,  thermal 
energy  is  transferred  from  the  electron  gas  to  the  molecules.  Even  as 
the  electron  density  approaches  zero,  75  percent  of  the  energy  released 
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during  recombination  is  carried  away  by  other  electrons  to  be  transferred 


to  the  gas  by  collisions.  Based  on  the  probe  I.D.  and  the  diameter  of 

-4 

the  four  bleed  holes,  the  residence  time  in  the  probe  is  about  10  sec- 

— _4 

onds.  Based  on  an  ion  flux  to  the  walls  of  1/4  n^  c^ , the  10  seconds 
is  plenty  of  time  for  all  the  ions  to  travel  to  the  wall  and  thermo- 
couple and  recombine,  releasing  the  energy  at  the  surface.  Based  on  the 
above  discussions,  it  is  reasonable  to  assume  complete  recombination 
within  the  total  temperature  probe. 

Another  mode  of  energy  loss  may  arise  from  existence  of  the 
metastable  states  of  Argon.  Following  the  method  of  Hinnov  and 
Hirschberg  [17],  the  number  density  of  the  metastable  states  will  be 
calculated.  The  assumption  is  that  the  metastable  states  are  in  transi- 
ent equilibrium  with  the  electron  temperature,  which  results  in  the 
equation: 
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Substituting  the  appropriate  values  based  on  the  experimental  data, 
we  have: 
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Thus  the  equilibrium  metastable  population  is  insignificant  when  compared 
to  the  electron  density,  and  the  contribution  to  the  energy  can  be 
ignored  in  the  calculations.  Calculations  of  the  power  radiated  from 
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the  gas  due  to  recombination  of  ions  and  electrons  using  experimental  and 
theoretical  recombination  coefficients  indicated  that  it  was  insignifi- 
cant compared  to  the  heat  input  to  the  gas. 

Since  an  R.F.  field  was  used  to  heat  the  gas,  magnetic  fields  were 
present  due  to  the  time  varying  electric  field.  It  then  is  necessary 
to  mention  the  effect  of  a magnetic  field  on  the  double  probe  operation. 
From  Chen  [18],  it  is  observed  that,  for  electrons  in  thermal  equilibrium 
and  an  ion  Larmor  radius  much  larger  than  the  probe  dimensions  (weak 
magneti  . field),  the  usual  theories  of  ion  collection  can  be  used.  In 
the  gas  flow  experiments,  the  electrons  appeared  to  be  Maxwellian,  the 
probe  was  far  enough  downstream  of  the  coil  so  that  the  magnetic  field 
was  weak,  and  ion  collection  was  the  probe's  mode  of  operation;  thus, 
under  the  test  conditions,  the  magnetic  field  had  no  noticeable  effect. 
R.F.  Energy  Addition 

The  energy  addition  in  the  experiments  was  provided  by  radio 
frequency  (R.F.)  induction  heating  using  a 5 KW  Lepel  R.F.  generator. 

The  electric  field  generated  by  the  R.F.  coil  initially  accelerates 
several  free  electrons  to  an  energy  where  they  can  ionize  Argon  atoms 
leading  to  more  electrons  and  an  avalanche  which  provide  an  ionized  gas 
or  plasma.  The  energy  transfer  for  heat  addition  can  be  considered  due 
to  electron  collisions  with  the  neutrals  and  ions  which  transfer  energy 
rather  than  ionizing  the  molecule  or  atom.  A simple  analysis  illustrates 
the  power  input  to  the  gas  by  electrons. 

The  momentum  equation  can  be  written  for  elastic  collisions  as: 


AFAPL-TR-76-57 


- e £«.  sin  out  - v me  d <jj 
<*+■  dt 


(112) 


The  collision  term  is  based  on  elementary  collision  theory  where  the 
electron  behaves  like  a very  light  elastic  sphere  which  loses  all  its 
momentum  in  any  given  direction  at  each  collision  with  a much  larger 
particle. 

Equation  (112)  can  be  written: 

C Oz  +■  V £>)  * = 


(113) 


Assuming  a particular  solution  of  equation  (113)  in  the  form 
ft*  - C cosuit  +Fsinwt  and  substituting  into  equation  (113)  solv- 
ing for  C and  F,  it  is  found  that: 
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the  drift  velocity  is  = dx/dt,  hence: 

V sin  tut 


V,  - efi 
d me 


- to  cos  git 
V2  +-to2 


in  cut 
f to2 


(115) 


the  current  density  is 

ne  clEo. 


J = 


me 


- q>  COS  tag 

V*  +■  u;2 


V Stn  tot 


•] 


V2  + *u2 

The  power  density  can  be  written  P 
applied  field.  Making  the  appropriate  substitutions,  it  is  found  that: 
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Multiplying  equation  (117)  by  dt,  dividing  by  the  period 
T = 1/f  = 2ir /.,  and  integrating  from  0 to  T,  the  average  power  per  unit 
volume  is  found  to  be: 


<P>  = 


_ nee2£* 


(118) 


2 me  v ^2 

2 

By  writing  the  kinetic  energy  of  an  electron  K.E.  = 1/2  mV^  , 
multiplying  by  dt,  dividing  by  the  period,  and  integrating  from  0 to  T, 
it  can  be  shown  that  the  average  kinetic  energy  of  an  electron  in  terms 


of  electron  volts  can  be  written: 
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Dividing  equation  (118)  by  equation  (119),  we  find: 
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or 
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The  electric  field  Ea  will  be  calculated  along  with  the  ion 

temperature  to  see  if  the  experimental  data  is  reasonable. 

<P>  is  determined  by  dividing  the  power  into  the  gas  by  the  volume 

2 2 

of  the  coil  region.  The  volume  is  given  by  Vol  = irL  (r-j  -r?  ) where  L 

is  the  coil  length  and  r^  and  r^  are  the  radius  boundaries  of  the  maximum 

temperature  region.  Substituting  the  appropriate  values,  we  find 
•6  3 

Vol  ~ 4 * 10  m . The  value  of  power  is  ' 100  watts,  which  is  about 

7 3 

the  maximum  in  the  experiments;  hence,  <P>  = 2.5  x 10  watts/meter  . 

The  term  <K.E.>  is  determined  by  the  measured  electron  temperature. 
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Using  the  maximum  values  from  the  1.6  area  ratio  experiments  for  electron 
temperature  (3  x 104  °K)  and  electron  density  (10^  particles/meter^) , we 
find  from  equation  (121)  that  v = 2.22  x 10  collisions  per  second  for 
= 2.5  x 107.  Substituting  the  above  values  back  into  equation  (119) 
and  solving  for  E , it  is  found  that  E,  2.3  volts/cm.  This  is  a rea- 
sonable  value  for  the  electric  field  in  the  plasma  based  on  data 
available  in  Brown  [19]. 

Next  let's  calculate  what  the  ion  or  neutral  temperature  will  be  in 


the  region  of  maximum  energy  input.  An  electron  in  an  elastic  collision 

will  transfer  about  2me/m^  of  its  energy  in  an  elastic  collision.  The 

-4 

residence  time  of  the  flow  in  the  coil  region  is  t 10  seconds.  The 
change  in  temperature  experienced  by  the  ions  and  neutrals  can  be  written; 

atl  = vt  rG  (122) 

Substituting  the  values,  it  is  found  that: 

^Tt ; = *26*/?  (123) 

Ti  - 5 SO  + 32  (.  = $7<S°/?  (124) 

550  is  initial  gas  temperature 


The  above  value  for  T.  compares  well  with  the  experimental  value  of  T. 

in  the  region  near  the  wall,  Ti  ~ 1000°R. 

w 

The  above  analysis  shows  that  the  measured  values  of  Tg,  ng,  and  T^ 
are  reasonably  correct  and  the  theory  can  describe  the  energy  addition 
process.  Using  Tg  to  calculate,  <K.E.>  is  reasonable  since 
ce  = "^8KTe/irme  = 8.7*10^  meters/second-,  and,  = 1/2  eEa/(mew),  which 
<K.E.>  depends  upon,  equals  7.4*1 0^  meters/second. 
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SECTION  IV 

EXPERIMENTAL  EQUIPMENT  AND  PROCEDURE 


Gas  Flow  Apparatus 

For  the  radio  frequency  (R.F.)  heating  of  the  inert  gas,  an 
electrically  non-conducting  material  was  needed  for  duct  fabrication. 
This  was  necessary  to  prevent  the  R.F.  water-cooled  copper  coil  (which 
surrounds  the  duct)  from  heating  the  duct  by  induction  currents  rather 
than  heating  the  gas.  Two  materials  were  used  to  make  the  ducts, 

Pyrex  glass  and  Boron  Nitride,  with  the  latter  being  used  in  the 
majority  of  experiments.  Boron  Nitride  is  a high  temperature,  easily- 
machined  refractory  material  with  a melting  point  of  2730°C;  and,  it  is 
a good  electrical  insulator.  Pyrex  glass  ducts  were  abandoned  after 
the  first  few  runs  because  they  became  plastic  during  high  energy 
inputs  to  the  gas.  The  gas  chosen  for  the  experiments  was  Argon,  since 
it  is  easily  ionized  and  relatively  inexpensive.  Ionization  of  the  gas 
is  necessary  so  that  the  electric  fields  generated  by  the  R.F.  coil  can 
induce  electric  currents  which  result  in  gas  heating. 

The  power  required  for  heating  the  gas  was  provided  by  a 0.5 
Kilowatt  (KW)  11  megacycle  R.F.  power  supply  or  a LEPEL  5 KW  4-8 
megacycle  power  supply;  the  latter  was  used  in  the  majority  of  the 
experiments.  The  R.F.  power  supply  was  connected  to  the  R.F.  coil  as 
shown  in  Figure  (20).  A Jennings  15  KV  0-25  uuF  variable  capacitor 
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was  used  to  tune  the  R.F.  duct  coil  to  the  power  supply.  This  tuning 
permitted  one  to  find  the  breakdown  point  (ionization  occurs)  and  to 
adjust  for  the  most  efficient  power  transfer.  The  copper  coil  wrapped 
around  the  duct  was  10-11  turns  of  1/8"  O.D.  water-cooled  tubing  of 
about  one  microhenry  inductance. 

To  ionize  and  heat  the  Argon  gas,  pressures  on  the  order  of  several 

millimeters  of  Hg  are  required.  To  achieve  these  low  pressures,  a 

vacuum  pumping  system  was  used  that  could  maintain  a back  pressure  of 

_2 

one  millimeter  of  Hg  at  an  Argon  flow  rate  of  10  lbm  per  second.  Some 
of  the  initial  experiments  used  a facility  with  about  half  the  above 
capability.  The  pumping  capacity  and  the  low  static  pressures  needed 
for  ionization  limited  the  total  pressure  in  the  nozzle  chamber  to  100 
mm  Hg;  this  total  pressure  restriction  limited  the  Reynolds  number  of 
the  flow.  To  achieve  the  high  Mach  number,  low  static  pressure  flow  in 
the  duct,  two  nozzles  were  fabricated  using  calculations  based  on  the 
analysis  of  Foelsh  [20].  His  analysis  provided  for  a parallel  flow  at 
the  nozzle  exit.  The  nozzles  had  design  exit  Mach  numbers  of  5 and  3.5 
with  an  exit  diameter  of  0.969  inches.  This  diameter  resulted  from  the 
fact  that  the  Pyrex  glass  tubing  used  initially  for  the  duct  had  an 
I.D.  of  0.969  inches.  The  Reynolds  number  based  on  the  duct  diameter 
for  the  Mach  3.5  nozzle  was  about  35,000  and,  as  mentioned  above,  this 
could  not  be  significantly  increased  due  to  ionization  and  pumping 
limitations. 

Figures  (21),  (22),  and  (23)  show  a schematic  and  pictures  of  the 
nozzle  and  duct  set-up.  Experiments  were  made  with  the  duct  directly 
connected  to  the  nozzle  and  separated  from  it,  as  shown  in  the 
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schematic.  To  allow  separations  of  the  nozzle  and  duct,  a small  chamber 
made  of  phenolic  was  fabricated  to  which  the  duct  and  nozzle  could  be 
attached.  The  small  chamber  had  two  ports  drilled  in  it  which  permitted 
a small  gas  bleed  into  one  port  and  pumping  on  the  other.  This  small 
chamber  allowed  for  a free  inlet  operation  of  the  duct  by  maintaining  a 
constant  pressure  in  between  the  nozzle  exit  and  duct  inlet.  The  pres- 
sure in  the  small  chamber  was  maintained  constant  by  having  a small  gas 
bleed  (about  four  percent  of  the  main  flow)  initially  into  the  chamber 
which  could  be  reduced  as  mass  spillage  occurred  at  the  duct  inlet  due 
to  heat  addition.  This  operation  permitted  maintenance  of  a constant 
Mach  number  in  front  of  the  duct  and  mass  spillage  due  to  heat  addition 
did  not  interfere  with  the  nozzle  flow.  Fischer  Porter  flow  tubes  were 
used  to  monitor  the  mass  flow  into  the  nozzle  and  the  small  connecting 
chamber.  Mass  flow  through  the  nozzle  based  on  the  measured  chamber 
total  temperature  and  pressure  agreed  within  two  percent  of  the  mass 
flow  indicated  by  the  flow  tube.  An  oil  bath  was  set  up  to  maintain  the 
temperature  of  the  Argon  entering  the  nozzle  at  a constant  90°F. 

Flow  Measurement  Probes 

The  Figures  (21),  (24),  and  (25)  illustrate  the  pressure  probes  and 
their  mounting  arrangement.  The  pressure  probes  were  constructed  from 
ceramic  tubing  of  0.98"  O.D.  and  0.05"  I.D.  Ceramic  instead  of  metal 
was  used  to  prevent  interaction  with  the  R.F.  coil.  The  probes  are 
about  5 inches  in  length  and  were  connected  to  stainless  steel  tubing 
which  goes  from  the  interior  of  the  vacuum  chamber  to  the  exterior  mount- 
ing system.  An  arrangement  was  made  so  that  two  probes  could  be  mounted, 
one  on  each  side  of  the  duct  and  nozzle  as  seen  in  the  above  figures.  A 
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threaded  rod  was  attached  to  the  stainless  steel  tubing  holding  the 
probe;  this  rod  passed  through  a threaded  knurled  nut  which  was  driven 
by  a low  rpm  motor  permitting  axial  motion  of  the  probe;  movement  was 
about  0.01"  per  second.  The  stainless  steel  tube  holding  the  probe 
passed  through  a double  0-ring  seal  in  the  flange  attached  to  the  vacuum 
chamber;  this  allowed  the  f?robe  to  be  rotated,  thus  permitting  radial 
measurements  to  be  made  iri  the  duct.  A securely-mounted  (though  still 
allowing  rotation)  slotted  stainless  tube  surrounded  the  one  to  which 
the  probe  was  attached;  this  acted  as  the  guide  for  a 1/4"  O.D.  tube 
mounted  90°  to  the  probe's  stainless  steel  tube.  A Tygon  tube  ran  from 
the  1/4"  O.D.  tube  to  a Datametrics  Barocel  Electrometer  Manometer,  which 
has  an  accuracy  of  _+  0.03°/  of  scale  reading.  Plots  of  static  and  total 
pressure  along  the  length  of  the  duct  for  a given  radial  position  were 
obtained  by  having  the  movable  probe  attached  to  a Bourns  0-5000  n 
Align-0-Plot  linear  resistor  and  a battery  which  was  connected  to  the  X 
axis  of  a X-Y  plotter.  The  0-10  volt  output  of  the  electronic  manometer 
was  connected  to  the  Y axis  of  the  X-Y  plotter.  See  Figure  (26)  for  a 
sketch  of  X-Y  plotter  hook-up.  Pressure  measurements , static  and  total, 
were  made  along  the  duct  at  radial  positions  of  O.OR,  0.19R,  0.38R,  0.57R, 
and  0.76R  from  the  center  line.  From  these  pressure  measurements , the 
Mach  number  along  the  duct  may  be  determined,  using  the  Rayleigh  pitot 
tube  formula. 

The  0.08"  O.D.  by  0.05"  I.D.  ceramic  tubes  were  used  for  the  probes 
because  Pyrex  tubing  became  plastic  at  the  higher  power  levels  when 
placed  near  the  duct's  wall.  One  total  probe  was  just  a flat-nosed 
ceramic  tube;  the  other  total  probe  had  a truncated  Boron  Nitride  5° 
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Figure  26.  Pressure  Measurement  Schematic. 
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half-angle  conical  tip  with  a 0.026"  hole.  The  static  probe  had  a 5° 
half-angle  1/2"  long  Boron  Nitride  conical  tip  with  four  holes  of  0.013" 
diameter  ninety  degrees  apart  at  a distance  of  1.0"  from  the  cone's  tip 
(the  holes  were  drilled  by  Bendix  Corp.  using  a cavitron  drill).  The 
total  temperature  probe  had  a Cromel -A1 umel  thermocouple  placed  in  the 
center  of  the  0.08"  O.D.  ceramic  tube  1/16"  from  the  inlet  with  four 
0.013"  diameter  holes  ninety  degrees  apart  1/8"  from  the  inlet  to  allow 
for  a small  mass  bleed  to  help  cancel  heat  conduction  losses.  The 
thermocouple  is  held  in  the  center  of  the  tube  by  a ceramic  spacer  to 
prevent  it  from  touching  the  tube  wall.  Figures  (24),  (25),  and  (27) 
show  the  probes  described  above.  Inclosing  the  thermocouple  provides  a 
low  velocity  flow  around  the  thermocouple  for  which  the  static  temperature 
is  essentially  the  total  temperature. 

For  use  in  conjunction  with  the  total  temperature  probe,  a Langmuir 
double  probe  was  made  to  determine  ionization  levels  so  that  corrections 
for  ionization  could  be  made  to  the  measured  total  temperature.  See 
Figures  (28)  and  (29)  for  illustrations  of  the  double  probe.  The  double 
probe  was  constructed  from  ceramic  tubing  and  0.005"  tungsten  wire.  The 
wires  were  spaced  about  1/3  inch  apart  and  were  about  1/4  inch  in  length. 
The  voltage  was  swept  across  the  wires  by  a variable  bi-polar  power  supply 
and  the  current  was  measured  across  a resistance  in  the  probe  circuit  by 
the  X-Y  plotter.  Reading  the  probe  voltage  on  the  other  plotter  axis 
resulted  in  a plot  of  probe  current  vs.  voltage,  see  Figure  (30)  for  a 
schematic  of  the  circuit.  The  static  and  total  pressures  were  measured 
first  for  all  experiments;  then,  the  pressure  probes  were  removed  and  the 
total  temperature  probe  and  double  probe  were  attached  to  the  two  probe 
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mounts.  The  flow  conditions  were  then  repeated  for  the  previous 
experiments  to  obtain  the  corresponding  total  temperatures  and  ioniza- 
tion levels  at  the  duct  exit. 

A collector  was  fabricated  and  placed  downstream  of  the  duct  inside 
the  vacuum  chamber,  see  Figures  (31)  and  (32).  The  collector  was  built 
to  allow  mixing  of  the  hot  and  cold  gas  of  the  jet  before  emerging 
through  four  holes  around  the  sides.  Thermocouples  were  placed  in  the 
flow  at  these  holes  to  measure  the  gas  temperature  before  and  after  heat 
addition.  The  collector  was  covered  with  aluminum  foil  to  reduce  the 
amount  of  radiation.  The  collector  temperature  was  monitored  by  thermo- 
couples at  the  locations  shown  in  Figure  (31).  This  collector  was  used 
in  conjunction  with  the  power  supply  meter  readings  and  mass  by-pass 
readings  to  assure  that  the  same  flow  conditions  were  repeated  when  the 
pressure  probes  were  removed  and  the  experiments  were  repeated  using  the 
total  temperature  and  Langmuir  double  probe.  The  collector  could  also 
be  used  for  a power  measurement  check  by  using  the  measured  temperature 
difference  before  and  after  heat  addition  and  the  mass  flow  rate.  A 
correction  to  the  power  input  could  be  made  for  collector  radiation 
loss  by  using  the  collector  temperature  and  its  surface  area  and  the 
emissivity  of  aluminum  (polished). 

Water  Table  Apparatus 

As  noted  in  the  introduction  and  discussed  earlier,  an  analogy 
between  water  flow  and  gas  flow  exists.  For  heat  addition  to  a gas  this 
implies  mass  addition  to  the  mainstream  water  flow.  To  accomplish  this 

simulation  the  set-up  illustrated  in  Figures  (33),  (34),  and  (35)  was  ' 

made  for  use  on  the  Air  Force  Institute  of  Technology's  water  table. 
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Figure  32.  Collector  As  Seen  Inside  The  Vacuum  Chamber. 

To  The  Right  Is  The  Port  In  Which  The  Collector 
Is  Placed  During  Gas  Flow  Tests. 
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The  water  table  is  a device  with  a glass  bottom  which  permits  water  flow 
from  one  end  to  the  other  at  various  adjustable  flow  rates.  The  glass 
bottom  enables  one  to  easily  observe  the  wave  patterns  on  the  water 
surface. 

A Mach  2 water  nozzle  was  fabricated  from  wood  based  on 
two-dimensional  wave  characteristics  for  a gas  having  a gamma  of  two. 

The  nozzle  was  placed  on  top  of  a 3/4"  thick  piece  of  plexiglas.  Down- 
stream of  the  nozzle  a plenum  was  cut  into  the  plexiglas  and  an  injection 
plate  with  numerous  holes  was  placed  over  this  plenum  flush  with  the 
surface.  Walls  for  the  simulated  duct  were  made  from  plexiglas  and  were 
attached  to  or  placed  a few  inches  downstream  of  the  nozzle  exit;  these 
walls  were  placed  on  both  sides  of  the  injection  plate  so  that  water 
could  be  injected  from  the  bottom  into  the  mainstream  flow  issuing  from 
the  water  nozzle.  This  injection  would  simulate  heat  addition  to  a gas 
flow. 

Photographs  were  taken  from  the  side  looking  at  the  duct  walls  for 
various  mass  injection  flow  rates.  These  photographs  permitted  observa- 
tion of  the  water  height  in  the  duct;  since  water  height  is  analogous 
to  the  gas  pressure,  the  water  height  permits  a qualitative  comparison 
with  the  gas  flow  tests.  These  water  table  and  gas  flow  experiments  are 
described  in  the  next  section. 
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SECTION  V 
EXPERIMENTS 

Water  Table  Experiments 

As  noted  earlier,  there  exists  a qualitative  analogy  between  gas 
flow  and  water  flow.  Therefore,  an  experiment  was  performed  on  the 
Air  Force  Institute  of  Technology  water  table  for  comparison  with  super- 
sonic gas  flow  energy  addition  tests.  These  experiments  qualitatively 
confirmed  the  phenomena  observed  in  the  gas  flow  experiments.  If  the 
nozzle  is  attached  to  the  duct,  heat  addition  or  mass  addition  causes  a 
pressure  rise  or  a water  level  rise  which  can  feed  upstream  in  the  sub- 
sonic portion  of  the  boundary  layer  into  the  nozzle,  thus  changing  the 
initial  conditions,  see  Figures  (36),  (37),  and  (38).  If  the  duct  is 
detached  from  the  nozzle,  the  pressure  rise  feeds  upstream  in  the  sub- 
sonic boundary  affecting  the  inlet  shock  and  causes  mass  spillage,  see 
Figures  (39),  (40),  and  (41).  The  effect  of  friction  is  visible  in 
Figures  (42)  and  (43)  where  the  length  of  the  duct  walls  are  increased 
(no  mass  addition  in  the  duct)  to  allow  for  a significant  friction 
increase.  The  effect  of  friction  is  similar  to  mass  or  heat  addition; 
a rise  in  pressure  or  water  level  occurs  which  can  be  transmitted 
upstream  through  the  subsonic  boundary  layer. 

For  a small  mass  flow  injection  there  was  a small  but  continuous 
rise  in  water  level  or  pressure  as  one  would  expect  in  an  ideal  weak 
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Figure  40 


Side  View  of  Free  Inlet  Duct 
Greater  Than  in  Figure  39. 


Figure  41 


Side  View  of  Free  Inlet  Duct 
Greater  Than  in  Figure  40. 
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Figure  44.  Static  Pressure  vs.  Axial  Position 
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detonation;  but,  with  increasing  mass  injection,  compression  and 
expansion  regions  occur  in  the  water  duct  as  was  observed  in  one  of 
the  initial  gas  flow  tests,  see  Figures  (38),  (41),  and  (44).  In 
most  cases  observed  in  gas  flow  experiments,  there  existed  a pressure 
rise  and  decrease  in  the  duct  as  in  Figures  (37)  and  (40).  In  all  the 
water  table  pictures  shown,  the  flow  remained  supersonic,  which  was 
verified  by  placing  a pencil  at  various  points  in  the  duct  and  observ- 
ing the  Mach  waves  generated  by  it. 

The  water  table  gave  a good  qualitative  picture  of  what  occurred  in 
the  gas  flow  energy  addition  experiments.  In  both  cases  the  pressure 
increase  along  the  duct  differed  from  the  gradual  continuous  increase 
predicted  by  the  ideal  theory.  As  will  be  shown  in  the  gas  flow  experi- 
ments this  difference  is  due  to  the  boundary  layer.  The  water  table 
facility  is  very  useful  in  practically  all  types  of  flow  problems.  If 
available,  it  should  be  used  to  gain  an  insight  into  a flow  problem 
before  moving  on  to  a more  expensive  and  time-consuming  gas  flow 
experiment. 

Gas  Flow  Experiments 

Three  ducts  were  used  for  the  supersonic  energy  addition  experiments 
having  area  ratios  of  1.0,  1.6,  and  2.0.  They  were  operated  at  three 
heat  input  levels.  See  Figure  (45)  for  a sketch  of  the  duct  with  an 
area  ratio  of  2.0.  Pressure  measurements  were  made  axially  and  radially 
along  the  ducts;  total  temperature  and  ionization  measurements  were  made 
radially  at  the  duct  exit.  The  temperature  and  ionization  measurements 
were  limited  to  the  duct  exit  because  an  insertion  of  about  1/2-inch 
into  the  duct  led  to  R.F.  interference  on  the  probe  output  rendering 
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useless  information;  in  addition,  the  diameter  of  the  double  probe  was 
too  large  for  insertion  without  seriously  disturbing  the  flow. 

Table  1 shows  the  results  of  total  temperature  measurements  made 
with  the  total  temperature  probe  for  several  nozzle  chamber  temperatures. 
No  energy  was  added  to  the  flow  in  the  duct.  The  nozzle  chamber  tempera- 
ture was  measured  on  the  center  line  as  were  the  duct  inlet  and  exit 
measurements.  As  the  chamber  total  temperature  was  increased,  the  probe 
measurement  error  decreased  to  about  2%  at  the  duct  exit. 

The  total  temperature  probe  had  a vent  hole  area  to  entrance  area 
ratio  (A  /A  ) of  0.27.  Goldstein  and  Scherrer  [21]  showed  that  the 
shielded  probe  recover  factor  r is  insensitive  to  Av/Ag  in  the  range  of 
0.2  to  0.8.  Also  in  the  range  of  M = 1.3  to  2.1  r was  a constant  value 
and  insensitive  to  angles  of  probe  inclination  up  to  10°.  Eber  [22] 
shows  a shielded  total  temperature  probe  with  bleed  holes  (similar  to  the 
one  used  in  this  effort)  that  is  insensitive  to  angles  of  attack  of  +10°; 
and,  in  the  range  of  M = 1 to  3.0,  the  recovery  factor  is  constant  at 
0.95.  Bontrager  [23]  determined  that  for  a ratio  of  entrance  length 
(distance  from  entrance  to  the  thermocouple  in  a shielded  probe)  to 
internal  diameter  equal  to  one  and  the  ratio  of  an  effective  friction 
length  to  entrance  length  greater  than  one,  then  the  measured  total 
temperature  approaches  closely  the  actual  value.  The  effective  friction 
length  is  that  length  equivalent  to  a 45°,  90°,  or  180°  turn  in  the  flow. 
Le/L  for  the  present  work  was  * 60.  Hence,  based  on  the  data  in  Table  1 
and  the  results  of  the  above  authors,  the  total  temperature  measured 
should  be  accurate  within  at  least  5%  of  the  actual  value.  However,  due 
to  the  recombination  of  ions  and  electrons  in  the  probe,  the  measured 
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value  must  be  corrected  for  recombination  energy  as  was  discussed  in  the 
double  probe  analysis  given  earlier.  The  error  due  to  this  recombination 
energy  correction  is  more  significant  than  the  error  in  the  total 
temperature  reading. 


TABLE  1 

TOTAL  TEMPERATURE  PROBE  MEASUREMENTS 


M,  = 3.1 

T°  - Measured  by  Probe 


Recovery  Recovery 


Nozzle  Chamber 
°R 

Duct  Inlet 
T°/T  l 

Factor 

r 

Duct  Exit 

t°/tJ. 

Factor 

r 

534.5 

0.982 

0.975 

0.974 

0.966 

550.5 

0.991 

0.984 

0.979 

0.974 

560.0 

0.995 

0.993 

0.982 

0.979 

580.0 

1.000 

1.000 

0.981 

0.977 

Initial  static  and  total  pressure  measurements,  using  the  probes 
discussed  earlier,  in  a diverging  duct  attached  to  a Mach  5 nozzle 
showed  large  errors  in  calculated  total  pressure  which  increased  with 
increasing  Mach  number  along  the  duct  axis.  Therefore,  a Mach  3.5  noz- 
zle was  fabricated  along  with  a constant  area  duct  with  three  pressure 
wall  taps  spaced  one  inch  apart.  The  purpose  of  this  duct  was  to 
observe  the  effect  of  the  pressure  probes  as  they  moved  down  the  duct. 
While  the  probe  was  moving  down  the  duct,  the  pressures  at  the  wall 
taps  were  recorded.  No  energy  was  added  to  the  duct  during  these  runs. 
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i Figure  (46)  shows  the  effect  of  the  conical  tipped  static  pressure  probe 

! at  the  wall  taps.  Figure  (47)  shows  the  effect  of  a flat  nosed  total 

probe  which  has  a considerable  upstream  influence  on  the  flow  due  to  the 
strong  shock  at  the  nose  interacting  with  the  boundary  layer.  The 
conical  tipped  static  probe  had  a small  upstream  influence.  Thus  it  was 
decided  to  fabricate  a total  pressure  probe  with  a truncated  conical  tip 
i to  reduce  the  upstream  influence  of  the  probe  which  introduces  errors 

| into  the  calculations. 

The  first  truncated  conical  tip  for  the  total  pressure  probe  was  0.3 
inches  long  and  had  a base  of  0.08"  diameter  through  which  a 0.013" 
diameter  hole  was  drilled  which  was  followed  by  a 0.025"  diameter  hole 
that  fed  into  the  0.05"  I.D.  of  the  ceramic  tube,  see  Figure  (24).  This 
probe  reduced  the  upstream  pressure  influence  significantly,  as  shown  in 
Figure  (48).  However,  when  the  measured  total  pressure  was  converted  to 
actual  total  pressure,  the  result  was  consistently  low,  whereas  the 
blunt-nosed  probe  gave  results  that  were  generally  high  when  compared  to 
the  chamber  total  pressure.  Checking  the  literature  on  total  probes, 
[24],  [25],  it  was  found  that,  according  to  the  presented  curves,  the 
dimensions  of  the  tip  hole  diameter  and  outside  diameter  of  the  probe  and 
the  Reynolds  number  were  such  that  the  probe  should  give  good  results. 

It  was  then  decided  to  double  the  hole  diameter  to  0.026"  all  the  way 
through  the  conical  tip,  keeping  the  outside  dimensions  the  same.  This 
correction  brought  the  calculated  total  pressure  much  closer  to  the 
nozzle  chamber  total  pressure.  The  maximum  errors  which  occur  in  regions 
of  rapid  (with  respect  to  axial  position)  pressure  variation  were  -20% 
and  +15%  compared  to  -50%  and  +200%  for  the  blunt-nosed  total  probe.  In 
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Figure  46.  Pressure  for  Wall  Taps  & Static  Pressure  Probe. 
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regions  of  slow  (with  respect  to  axial  position)  pressure  variation,  the 
calculated  total  pressure  ranged  from  -6"  to  +6%  of  the  nozzle  chamber 
total  pressure. 

To  illustrate  the  difficulty  of  determining  the  actual  total 
pressure  from  measured  total  and  static  pressure,  we  can  write  the 
differential  equation: 


ACT  _ y 

y- 


flCT 


Jll 

P° 


clP 

i? 


(125) 


where  P & PQ  is  the  measured  static  and  total  pressure  behind  the  probe 
shock  wave.  Since  y/y-1  = 2.5,  then  the  combination  of  the  errors  in 
measured  total  and  static  pressure  can  cause  significant  errors  in 
calculated  total  pressure.  For  example,  a 5%  error  in  P°  and  a -5,  error 
in  P will  result  in  an  error  of  20%  in  P0^-]--  Comparing  the  static  pres- 
sure as  measured  by  the  probe  with  that  of  the  wall  taps,  it  was  observed 
that  probe  and  wall  tap  agreed  within  a percent  when  in  a region  of 
relatively  slow  pressure  variation  and  the  holes  of  the  probe  were  in 
line  with  the  wall  tap. 

Using  the  probes  and  three  ducts  discussed  above,  measurements  were 
made  for  zero  power  and  three  power  settings  for  each  duct.  Pressure 
measurements  were  made  axially  along  the  duct  at  the  center  line,  0.19R, 
0.38R,  0.57R,  and  0.76R  from  the  center  line.  These  pressures  were 
recorded  on  an  X-Y  plotter  continuously  and  a digital  printout  every  1/4 
inch.  The  total  temperature  and  ionization  measurements  were  made  at 
11  positions  radially  from  the  center  line  at  the  duct  exit.  These 
measurements  were  also  recorded  on  an  X-Y  plotter.  The  probes  moved  at 
about  0.01  inches  per  second;  but,  even  at  this  slow  speed  a few  seconds' 
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pause  was  necessary  at  each  1/4"  interval  to  allow  the  static  pressure 
to  reach  equilibrium.  Figure  (49)  shows  the  probes  in  the  jet  downstream 
of  the  duct  during  a power-on  run. 

A computer  program  was  utilized  to  analyze  the  data  and  to  plot 
properties  vs.  radius  or  average  properties  vs.  axial  position.  Equa- 
tions (polynomials)  as  a function  of  radius  were  developed  to  satisfy 
the  data  points  and  the  boundary  conditions  at  the  duct's  center  and 
wall.  These  equations  listed  as  subroutines  in  the  computer  program 
could  then  be  utilized  to  determine  averages  by  integration  and  be 
called  by  the  plot  routine  to  develop  the  radial  or  axial  plots.  Figures 
(50)  through  (64)  show  the  experimental  results  for  the  2°  24'  half-angle 
diverging  duct.  Tables  2 through  10  compare  the  calculated  average 
properties  of  the  data  for  all  three  ducts  with  the  analytical  program 
based  on  average  properties,  which  was  discussed  in  Section  I.  The  input 
for  the  theoretical  program  was  the  average  initial  and  final  Mach  num- 
bers from  the  experimental  data  taken  at  inlet  and  exit  of  the  duct;  then 
other  parameters  were  calculated  for  comparison.  Average  properties  for 
the  data  were  calculated  two  ways;  one  was  integrated  by  RdR  and  divid- 

9 

ed  by  R“/2,  or,  if  a parameter  is  a function  of  two  or  more  variables, 
an  average  value  was  obtained  by  using  the  average  value  of  each 
variable  involved;  for  example,  p = <P>/R<T>.  For  the  temperature,  an 
average  by  TmRdT  divided  by  mRdT  was  also  made. 
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NOTES  FOR  FIGURES  (50)  THROUGH  (64) 


Curves  are  for  a duct  with  area  ratio  of  1.6. 

Total  pressure  and  total  temperature  ratios  are  referenced  to  the 
nozzle  chamber  conditions.  Other  parameter  ratios  are  referenced  to 
the  value  at  the  duct  inlet  center  line  before  energy  addition. 


a ->  degree  of  ionization 

(a  = 0 unless  noted  otherwise) 

I -*■  refers  to  position  along  the  duct  inlet  (1)  to  exit  (20)  in 
intervals  of  1/4" 

N refers  to  power  setting  or  energy  input  given  by  total 
temperature  ratio 
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FIGURE  55  . TOTAL  PRESS.  RATIO  VS. 
RADIAL  POSITION 


POI-  100  flrt.  MG. 
TO  - 550.0  R 
PI  = 2.40  MM.  MO. 
T 1 = 130. B R 
Ml  = 3. 10 


A VE  TOTAL  PRES  RATIO  <P0/P01> 

.27  0.35  0.43  0.51  0.59  0.67 


AF  vPL-  ' 


N 1 2 3 4 POt=  100  MM.  HO. 

1.0  1.03  1.1?  1.19  T02/T01  TO  = 560.0  R 

PI  = 2.40  MM.  HO. 

0.0  0.16  0.59  3.57  X MASS  8Y  PASS 

T 1 = 130.8  R 

Ml  = 3. 10 

I 2C 


. c.  u u.tu  u . u u 

AXIAL  POSITION  Z/L 
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To  determine  the  friction  factor  for  the  constant  area  duct,  the 
experimental  average  inlet  and  exit  Mach  numbers  for  cold  flow  were 


used  in  the  ideal  constant  area  friction  factor  equation: 


j±L  _ ± 

P 


(126) 


This  equation  gave  a value  of  f = 0.0003.  For  the  diverging  duct  the 
friction  factor  term  f/tan9  was  determined  by  using  the  initial  and 
final  average  Mach  numbers  determined  from  the  experiment  for  cold  flow 
and  varying  f/tano  until  the  proper  area  ratio  was  obtained.  The  value 
for  f/tan8  was  = 0.093;  however,  to  establish  reasonable  agreement 
between  the  measured  and  calculated  total  temperature  ratios  and  the 
corresponding  area  ratio,  it  was  necessary  to  correct  the  value  of 
f/tane  by  multiplying  it  by  the  square  root  of  the  ratio  of  the  wall 
temperature  after  heat  addition  to  that  before  heat  addition.  The 
reason  for  the  above  is  that  in  laminar  flow,  which  was  the  experimental 
condition,  the  friction  factor  f is  proportional  to  1/Re  which  makes  it 
proportional  to  the  viscosity  y,  f is  the  wall  friction  term.  For  a 
monatomic  gas  such  as  Argon  the  viscosity  is  proportional  to  TT;  thus , 
as  the  gas  temperature  increases  at  the  wall,  one  can  expect  the  wall 
friction  term  to  increase  by  TT.  The  values  marked  by  an  asterisk  in 
Tables  2 through  10  represent  those  values  close  to  the  corrected 
f/tane  term. 

The  values  of  N and  n for  use  in  the  theoretical  analysis  were 
determined  from  the  experimental  data  as  follows.  Total  temperature 
ratio  values  were  determined  for  the  following  cases: 


■ 
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,L 


1 • 


i 

I 


i 


1 . ^VRdR/jf^RdR 

2.  (^^T°RdR/^^RdR) . (i  -►  refers  to  ionization  correction) 

3.  ^RT°mRdR/^RmRdR 

4.  (jfRT°mRdR/^RmRdR). 

5.  2 r^rOw/A  RdR  •*  converted  to  total  temperature 

(Pow/A  ->-  power  per  unit  area) 

6.  (2* /£RPow/A  RdR)^  > converted  to  total  temperature 

Power  per  unit  area  written  in  terms  of  average  properties  and  converted 
to  total  temperature  where  averages  are  determined  as  follows: 

P -*■  parameter 

7.  ^RPRdR/^RRdR 

8.  (-jfRPRdR/>jrRRdR)i 

9.  (J^RPmRdR//'RmRdR) 

10.  (^RPmRdR/jTRmRdR). 
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Knowing  the  area  ratio  and  using  the  above  total  temperatures  ratios,  N 
was  calculated  from  T2/T°  = (A2/A-|)^.  From  the  experimental  flow  rate 
measurements,  n was  determined  by  m^/n^  = (A2/A-|  )n.  All  these  cases 
are  listed  in  Tables  2 through  10.  For  the  constant  area  case,  n and 
N were  determined  by  letting  A9/A^  be  a number  only  slightly  greater 
than  one.  Table  11  shows  ionization  values  and  electron  temperatures 
for  the  constant  area  duct  experiments  as  a function  of  radius. 
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TABLE  11 

CONSTANT  AREA  DUCT 


Double  Probe  Data  for  Three  Total  Temperature 
Ratios  as  Function  of  Radial  Position 


Ionization 


X 104 


T°/T°  T°/T°  T°/T°  T°/T°  T°/T° 
1.05  1.10  1.18  1.05  1.10 


0.0  0.092  0.531  0.665  2.58  2.44  3.87  0.023  0.121  0.161 


0.095  0.115  0.706  0.811  2.00  1.800  3.50  0.029  0.165  0.208 


0.1900.1801.2351.0781.37  1.12  3.06  0.0430.2780.271 


1.00  1.13  3.95  0.297  0.601  0.376 


0.805  1.10  2.10  0.398  1 .079  1.236 


1.18  1.18  1.77  0.327  1 .752  2.047 


1 .47  1 .84  0.738  2.305  3.141 


1 .53  1.90  0.911  2.521  3.669 


1.66  2.04  1.026  2.234  3.417 


1.79  2.19  1.177  1.989  3.183 


2.472  2.995 


0.0 


1.11 

1 .20 

2.10 

0.0 

0.0 

0.0 
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SECTION  VI 

SUMMARY  AND  CONCLUSIONS 

Experiments 

The  flow  analyzed  was  heat  addition  to  a supersonic  flow  of  Argon 
gas  passing  through  a 4 3/4-inch  long  Boron  Nitride  (B.N.)  duct  with  a 
one-inch  diameter  inlet.  The  duct  was  located  one-half  inch  downstream 
of  an  axisymmetric  nozzle.  The  nozzle,  designed  for  parallel  flow  at 
the  exit,  had  a one-inch  exit  diameter  and  a stagnation  chamber  pres- 
sure and  temperature  of  100  mm  Hg  and  550°R.  The  center  line  Mach 
number  at  the  duct  entrance  was  approximately  3.1,  while  the  average 
across  the  duct  diameter  was  about  2.5.  The  flow  inside  the  duct  con- 
tained compression  and  expansion  regions  which  were  a result  of  the 
flow  not  entering  exactly  parallel  to  the  duct  inlet  lip.  These  com- 
pression and  expansion  regions  remained  essentially  the  same  during 
heat  addition  to  the  gas.  The  static  pressure  in  the  gap  between  the 
nozzle  exit  and  duct  entrance  was  maintained  constant  by  a small  mass 
flow  bleed  which  could  be  decreased  to  compensate  for  any  mass  flow 
bypass  caused  by  heat  addition  to  the  gas  flow. 

Heat  was  added  to  the  supersonic  duct  flow  by  means  of  a one  and 
one-half  inch  long  R.F.  copper  coil  wrapped  around  the  B.N.  duct  and 
connected  to  a five  Kilowatt  Lepel  R.F.  power  supply.  Experiments 
were  made  with  ducts  having  area  ratios  of  1.0,  1.6,  and  2.0.  For 
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each  duct,  three  heat  addition  values  were  investigated.  The  necessary 
measurements  were  made  with  a static  pressure  probe,  a total  pressure 
probe,  a total  temperature  probe,  and  a Langmuir  double  probe.  The 
probes  were  made  of  ceramic  to  avoid  interference  with  the  R.F.  dis- 
charge and  to  withstand  the  high  temperature  environment.  Static  and 
total  pressure  measurements  were  made  at  five  different  radial  posi- 
tions along  the  length  of  the  duct.  The  double  probe  and  total  tempera- 
ture measurements  were  made  only  at  the  duct  exit  because  the  double 
probe  was  too  large  to  insert  without  drastically  changing  the  flow  and 
the  thermocouple  began  to  interact  with  the  R.F.  coil  when  inserted 
about  one-half  inch  into  the  duct  exit.  The  double  probe  was  used  to 
measure  ionization  levels  so  that  a correction  for  ion-electron  recom- 
bination could  be  made  to  the  total  temperature  probe  measurement. 

A water  table  experiment  was  performed  to  determine  if  the  hydraulic 
analogy  would  qualitatively  predict  the  flow  phenomena  for  heat  addition 
to  a supersonic  gas.  Heat  addition  was  simulated  by  water  mass  flow 
addition  to  a stream  of  supersonic  water.  The  mass  addition  was  accom- 
plished by  injection  through  a porous  plate  beneath  the  main  water 
stream.  A properly-contoured  wood  nozzle  provided  the  Mach  2 water  flow 
and  the  duct  walls  were  made  of  plexiglas  for  visual  observation  of  the 
water  height.  The  results  of  this  experiment  qual i tati vely  showed  the 
same  pressure  (water  heiqht)  profiles  along  the  duct  as  observed  in 
the  gas  experiments. 

Radial  and  Axial  Variations  for  Gas  Experiments 

The  graphs  presented  in  Figures  (50)  to  (64)  illustrate  the  radial 
and  axial  property  variations  for  the  1.6  area  ratio  duct.  The  Mach 
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number  profiles  are  fairly  flat  out  to  about  0.7  of  the  radius  and  show 
the  proper  trend  for  heat  addition,  that  is,  decreasing  with  increasing 
heat  input.  The  graph  of  static  pressure  versus  axial  position  illus- 
trates two  important  points.  First,  the  pressure  rise  due  to  heat  addi- 
tion propagates  upstream  into  the  inlet  even  though  the  flow  is  super- 
sonic and  second,  the  pressure  first  increases  and  then  decreases 
within  the  duct  for  the  two  highest  heating  values.  The  upstream  pres- 
sure rise  can  only  be  due  to  the  subsonic  portion  of  the  boundary  layer 
which  permits  the  pressure  disturbance  from  heat  addition  to  propagate 
upstream  into  the  inlet  region,  causing  mass  flow  spillage.  The  pres- 
sure rise  and  decrease  within  the  duct  can  be  explained  if  a separation 
region  occurs  within  the  duct  and  reattaches  near  the  duct  exit.  The 
graphs  of  density  ratio  and  mass  flow  ratio  in  Figures  (62)  and  (64) 
show  that  density  and  mass  flow  both  decrease  near  the  wall  and  increase 
in  the  center  region  of  the  duct.  Therefore,  this  could  indicate  a 
converging  of  the  streamlines  caused  by  heat  addition  and  separation 
due  to  an  adverse  pressure  gradient.  From  Shapiro  [8],  Vol  II,  the 

pressure  ratio  necessary  for  separation  of  a laminar  boundary  layer 

2 

interacting  with  a shock  wave  can  be  written  P^/P]  = 1 + 1/2  yM^f 
(Reynolds  number).  For  the  conditions  in  the  experiment,  P^/P^  is  equal 
to  1 + .05  For  = 3.1,  we  find  P^/ P i 1.8  which  is  exceeded  in 

the  higher  heat  input  values  between  the  inlet  and  center  of  the  duct; 
thus,  separation  is  a possibility.  Separation  causes  an  effective 
decrease  in  flow  area  which  increases  the  static  pressure  in  a super- 
sonic flow.  A separation  thickness  of  about  0.065"  would  result  in  a 
pressure  ratio  of  about  two  as  a result  of  the  decrease  in  flow  area. 
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The  magnetic  pressure  force  was  estimated  to  be  on  the  order  of  0.001 
mm  Hg  near  the  wall  and  hence  should  not  be  a significant  factor  in 
flow  separation  since  the  duct  pressure  is  on  the  order  of  4 mm  Hg  or 
higher. 

The  wave  pattern  observed  in  the  duct  could  be  the  result  of  the 
interaction  of  the  oblique  shock  wave  with  a laminar  boundary  layer  as 
discussed  in  Shapiro,  page  1141,  Vol . II  [8].  The  interaction,  how- 
ever, implies  a separated  zone  which  occurs  at  a pressure  ratio  less 
than  that  mentioned  above.  According  to  Young  [26],  a pressure  rise  due 
to  a shock  wave  impinging  on  a laminar  boundary  layer  can  extend  100 
boundary  layer  thicknesses  upstream  of  the  intersection  point.  This 
upstream  influence  will  occur  whether  the  boundary  layer  separates  or 
not.  Thus,  the  pressure  disturbance  from  heat  addition  can  likewise 
propagate  upstream  a substantial  distance  due  to  the  subsonic  part  of 
the  boundary  layer.  This  upstream  pressure  propagation  causes  a thicken- 
ing of  the  boundary  layer  which  effectively  increases  the  flow  area  and 
hence  increases  the  duct  static  pressure.  The  large  pressure  rise  and 
decrease  within  the  duct  during  heat  input  could  be  due  to  a combina- 
tion of  two  phenomena.  The  pressure  variation  can  be  caused  by  the  flow 
separation  discussed  previously;  and,  secondly,  the  heating  and  expan- 
sion of  the  gas  in  the  R.F.  coil  region  forces  mass  towards  the  center 
of  the  duct  resulting  in  an  effective  area  decrease.  Following  this 
rapid  heating,  radial  conduction  of  heat  occurs,  thus  relieving  the 
pressure. 

The  figures  for  static  and  total  temperature  show  that  most  of  the 
temperature  rise  occurs  within  a region  out  to  about  1/2  the  duct 
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radius  from  the  wall.  This  region  occurs  because  the  R.F.  discharge 
tends  to  shield  the  interior  flow  from  the  electric  and  magnetic  fields 
of  the  coil.  This  shielding  can  be  considered  similar  to  the  skin 
effect  in  a metal  conductor.  Figure  (63)  for  power  vs.  radial  position 
illustrates  that  even  though  the  center  region  of  the  duct  receives 
little  R.F.  heating,  the  power  per  unit  area  increases  largely  due  to 
the  heated  gas  near  the  wall  forcing  mass  flow  into  the  center  region 
of  the  flow. 

From  the  plot  of  average  total  pressure  along  the  duct  axis,  it  is 
observed  that  the  total  pressure  tends  to  decrease  along  the  duct  with 
heat  addition  and  then  increases  again  before  reaching  the  exit.  This 
trend  is  due  to  errors  in  the  measured  total  and  static  temperatures. 

As  mentioned  in  the  section  entitled  Gas  Flow  Experiments,  a -51  error 
in  measured  static  pressure  and  a 5%  error  in  impact  pressure  can  lead 
to  a 20»  error  in  calculated  total  pressure.  Thus  the  fluctuations  of 
average  total  pressure  along  the  axis  is  not  surprising.  Figure  (55) 
shows  total  pressure  as  a function  of  radius  at  the  duct  exit;  the 
total  pressure  decreases  with  increasing  heat  input  as  it  should  in  a 
supersonic  flow  with  heat  addition. 

Average  Properties  for  Gas  Experiments 

A quasi-one-dimensional  theory  based  on  average  properties 
(averaged  with  respect  to  radius)  was  developed  using  the  assumptions 
( T^  / T°  ) = (A^/A^)^;  m^/m-i  = (A^/A^)0,  where  N and  n are  variable 
parameters.  The  above  two  assumptions  are  reasonable  since  eight 
special  flow  conditions  can  be  described  by  them  (see  section  entitled 
Flow  Equations).  In  general,  one  can,  by  using  N and  n as  variable 
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parameters  (and  including  wall  friction  by  means  of  the  friction  factor, 
f),  analyze  flows  in  axisymmetric  ducts  with  area  change,  energy 
exchange,  and  variable  mass  flow  occurring  simultaneously.  The  theory 
was  developed  to  compare  with  the  experimental  results. 

Tables  2 through  10  show  the  calculated  average  experimental  results 

with  those  calculated  by  the  quasi-one-dimensional  theory.  The  N and  n 

required  were  determined  from  the  total  temperature  averages  discussed 

in  the  section  entitled  Gas  Flow  Experiments  and  the  measured  mass 

flows.  The  initial  and  final  Mach  numbers  needed  for  the  theoretical 

inputs  were  based  on  the  average  inlet  Mach  number  before  heat  addition 

and  the  average  Mach  number  at  the  exit  after  heat  addition. 

fR  o*  fR  . 

It  is  seen  from  the  above  tables  that  averaging  by  J ^ T mRdR/_/Q  mRdR 
was  more  appropriate  than  averaging  by  J^RT°RdR/jTRRdR  as  assumed  in  the 
quasi-one-dimensional  theory.  The  reason  is  that  an  energy  change 
(which  is  related  to  the  total  temperature  or  static  temperature)  is  a 
unit  mass  property  and  when  significant  radial  variations  exist,  a sub- 
stantial error  can  be  obtained  by  averaging  by  RdR.  The  static  tempera- 
ture also  agreed  better  with  theory  when  averaged  by  mRdR.  Total  tem- 
perature measured  by  the  probe  can  sense  the  energy  released  by  electron- 
ion  recombination;  this  energy  release  was  taken  into  account  by  the 
ionization  level  determined  from  the  Langmuir  double  probe;  and, 
averages  with  and  without  the  ionization  correction  are  shown  in 
Tables  2 through  10.  In  most  cases,  correcting  for  ionization  gave 
closer  agreement  with  the  theoretical  values. 


The  experimental  static  temperature  was  always  higher  than  the 
theoretical  value  which  indicates  the  total  temperature  probe  probably 
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read  too  high  a value;  this  difference  in  temperature  may  be  due  to 
the  error  in  the  ionization  level  correction,  which  could  be  off  by  an 
order  of  magnitude  of  two.  The  static  pressures  agreed  fairly  well 
between  theory  and  experiment  and  surprisingly  the  total  pressure  also 
agreed  well  in  many  cases.  The  value  of  f/tane  for  the-  diverging  area 
ducts  was  ' 0.093  as  determined  from  the  cold  flow  (no  heat  addition) 
condition.  The  values  of  Tables  5 through  10  were  obtained  by  varying 
f/tane  until  the  correct  area  ratio  was  obtained  for  a given  duct  and 
heat  addition  value.  The  correct  area  ratio  was  obtained  in  the 
majority  of  cases  for  f/tane  close  to  0.093  multiplied  by  the  tempera- 
ture correction  discussed  in  the  section  entitled  Gas  Flow  Experiments 
and  corresponded  to  averages  based  on  mRdR  and  corrected  for 
ionization. 

Figure  (65)  illustrates  the  average  total  pressure  variation  with 
increasing  average  T°/T°  and  increasing  area;  the  trend  is  correct, 
decreasing  for  increasing  T^/T®  and  decreasing  for  increasing  area 
which  implies  energy  addition  at  a higher  Mach  number.  Figure  (66) 
illustrates  the  mass  flow  ratio  bypass  with  heat  addition  for  constant 
area  and  diverging  ducts;  again  the  trend  is  correct,  the  mass  bypass 
is  decreased  with  increasing  area  for  a given  energy  input  because  of 
the  pressure  relieve  offered  by  a larger  area  duct. 

Conclusions 

A steady,  radially  non-uniform  weak  detonation  wave  corresponding 
to  point  c of  Figure  (1)  was  simulated  by  R.F.  heating  of  a supersonic 
Argon  gas  flow.  The  quasi-one-dimensional  theory  based  on  average 
properties,  where  T^/T®  = (A^/A^)^  and  m^/m^  = (A^/A^ )n  can 
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Figure  65.  Total  Pressure  Ratio  vs.  Total  Temperature  Ratio 
for  Constant  Area  and  Diverging  Ducts. 
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Figure  66.  Mass  Flow  Ratio  Caused  by  Heat  Addition  vs. 

Total  Temperature  Ratio  for  Constant  Area 
and  Diverging  Ducts. 
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satisfactorily  describe  the  weak  detonation  mode  for  constant  and 
varying  area  ducts,  even  though  the  flow  before  and  after  heat  addition 
is  not  one-dimensional  and  the  energy  inputs  are  non-uniform  across  the 
radius.  The  theory  works  when  the  initial  average  Mach  number  <M^ > is 
based  on  the  initial  inlet  condition  before  heat  addition.  When  deter- 
mining an  average  temperature  across  the  radius  for  comparison  with  the 
theory,  the  best  method  is  by  ^TmRdR  /JjB  RdR.  For  small  variations  in 
the  radial  direction,  averaging  by  RdR  is  appropriate. 

The  discrepancies  that  appear  between  some  of  the  experimental  and 
calculated  properties  are  due  to  measurement  errors  which,  however 
slight,  can  result  in  significant  errors  for  quantities  that  are  computed 
from  measured  properties;  this  was  shown  earlier  for  the  case  of  total 
pressure.  More  accurate  measurements  should  result  in  better  agreement 
between  theory  and  experiment.  Energy  addition  to  a supersonic  flow  in  a 
real  situation  can  cause  a pressure  rise  upstream  of  the  energy  input 
region  due  to  the  subsonic  portion  of  the  boundary  layer.  The  extent  of 
this  influence  depends  on  the  Reynolds  number;  and,  if  it  is  small  enough, 
significant  mass  bypass  can  be  expected  at  the  inlet  as  was  shown  in  the 
experiments.  The  average  Mach  number  in  the  duct  approached  1.0  only  in 
one  instance,  and  in  all  other  cases  was  above  1.0;  hence,  choking  of  the 
flow  was  not  responsible  for  the  upstream  influence  which  was  observed. 

The  water  table  experiments  qualitatively  exhibited  the  same  phenomena 
seen  in  the  gas  flow  tests,  thus  validating  its  use  for  gaining  an 
insight  into  gas  flow  energy  addition.  An  electric  double  probe  was 
built  to  me  sure  ionization  levels  in  a supersonic  gas  stream  and  it 
functioned  well.  The  equation  developed  for  the  degree  of  ionization 
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was  derived  by  analyzing  the  theories  of  Cohen,  Waymouth,  and  Lam  for 
dense  plasmas.  It  was  shown  that  all  three  theories,  although  of  dif- 
ferent initial  form,  reduced  to  the  same  expression  for  the  case  of  ion 
collection  at  large  negative  potentials  and  electron  temperatures  much 
greater  than  ion  or  neutral  temperatures. 

When  a practical  supersonic  combustion  engine  is  developed,  it  may 
rely  on  the  combustion  or  energy  addition  occurring  near  the  wall  region 
as  in  these  experiments;  so,  the  test  results  obtained  here  should  pro- 
vide an  insight  to  the  problems  to  be  encountered  in  high  altitude 
supersonic  combustion.  As  was  mentioned  in  the  Introduction,  supersonic 
energy  addition  can  be  by  a subsonic  diffusion  flame  propagating  in  a 
supersonic  stream  or  by  a weak  detonation  wave  which  travels  at  a super- 
sonic speed  relative  to  both  the  heated  and  unheated  gas.  The  latter 
is  difficult  to  obtain  with  a chemical  reaction,  but  was  simulated  in 
these  experiments  with  R.F.  heating.  However,  energy  addition  by  R.F. 
was  not  assumed  to  be  a practical  means  of  energy  transfer,  but  only  a 
method  to  achieve  a weak  detonation  without  resorting  to  a reacting  gas 
mixture.  To  be  practical,  R.F.  energy  addition  efficiency  must  be 
greatly  increased  over  the  10%  which  was  measured  in  a "dummy"  load  test. 

The  effort  has  added  to  the  knowledge  of  energy  exchange  with  a 
flowing  gas  by  analyzing  non-uniform  (radially),  low  Reynolds  number 
(35,000),  varying  area  supersonic  energy  addition.  The  major  energy 
input  was  restricted  to  a region  near  the  wall  of  the  duct.  A future 
effort  may  be  one  that  adds  energy  to  the  center  of  the  duct  by  passing 
an  electric  arc  along  the  center  line.  This  would  add  the  energy  in  a 
region  of  high  Mach  number  compared  to  the  wall  region;  for  the  same 
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energy  input,  it  would  be  interesting  to  see  if  the  results  are  the  same 
as  in  the  present  work.  If  the  results  are  equivalent,  then  this  would 
help  confirm  the  averaging  method  proposed  and  the  quasi -one-dimensional 
theory  outlined  earlier.  Another  area  that  could  be  explored  is  a 
non-steady  theory  that  can  predict  for  supersonic  heat  addition,  the 
upstream  pressure  rise  which  propagates  in  the  subsonic  portion  of  the 
boundary  layer;  in  the  ideal  supersonic  theory  there  is  no  upstream 
influence  as  long  as  no  choking  occurs.  The  achievement  of  a steady 
weak  detonation  wave  in  a chemically-reacting  gas  mixture  must  wait  the 
development  of  a highly  efficient,  low  power  mechanism  (perhaps  laser 
radiation)  which  will  initiate  and  sustain  the  chemical  reaction  in  a 
small  region. 
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